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We calculate the explicit equations of motion for non-spinning compact objects to 2.5 post- 
Newtonian order, or 0(v/c) 5 beyond Newtonian gravity, in a general class of scalar-tensor theories 
of gravity. We use the formalism of the Direct Integration of the Relaxed Einstein Equations 
£T") ■ (DIRE) , adapted to scalar-tensor theory, coupled with an approach pioneered by Eardley for incor- 

porating the internal gravity of compact, self-gravitating bodies. For the conservative part of the 
motion, we obtain the two-body Lagrangian and conserved energy and momentum through second 
post-Newtonian order. We find the 1.5 post-Newtonian and 2.5 post-Newtonian contributions to 
( — i ■ gravitational radiation reaction, the former corresponding to the effects of dipole gravitational ra- 

diation, and verify that the resulting energy loss agrees with earlier calculations of the energy flux. 
For binary black holes we show that the motion through 2.5 post-Newtonian order is observationally 
identical to that predicted by general relativity. For mixed black-hole neutron-star binary systems, 
the motion is identical to that in general relativity through the first post-Newtonian order, but 
deviates from general relativity beginning at 1.5 post-Newtonian order, in part through the onset 
of dipole gravitational radiation. But through 2.5 post-Newtonian order, those deviations in the 
motion of a mixed system are governed by a single parameter dependent only upon the coupling 
a* constant ujq and the structure of the neutron star, and are formally the same for a general class of 

5^ ■ scalar-tensor theories as they are for pure Brans-Dicke theory. 
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INTRODUCTION AND SUMMARY 



The anticipated detection of gravitational waves by a network of ground-based laser-interferometric observatories 
promises a new way of "listening" to the universe in the high-frequency band. A future space-borne interferometer 
would open the low-frequency band and pulsar timing arrays may soon begin exploring the nano-Hertz region of the 
gravitational-wave spectrum. In addition to providing a wealth of astrophysical information, these observations also 
hold the promise of providing tests of Einstein's theory of general relativity in the strong-field, dynamical regime. 
. — i ■ The "inspiralling compact binary" - a binary system of neutron stars or black holes (or one of each) in the late 
stages of inspiral and coalescence - is a leading potential source for detection. Given the expected sensitivity of the 
ground-based interferometers, stellar-mass compact binaries could be detected out to hundreds of megaparsecs, while 
for a space interferometer, inspirals involving supermassive black holes could be heard to cosmological distances. 

In order to maximize the detection capability and the science return of these observatories, extremely accurate, 
theoretically generated "templates" for the gravitational waveform emitted during the inspiral phase must be available. 
This means that correction terms in the equations of motion and gravitational- wave signal must be calculated to high 
orders in the post-Newtonian (PN) approximation to general relativity, which, roughly speaking, is an expansion in 
powers of v/c ~ (Gm/rc 2 ) 1 ^ 2 (for a review and references see [1]). Contributions to the waveform from the merger 
phase of the two objects and from the "ringdown" phase of the final vibrating black hole also play an important role. 

The detected gravitational- wave signals can also be used to test Einstein's theory in the radiative regime, particularly 
for waves emitted by sources characterized by strong-field gravity, such as inspiraling compact binaries. One way to 
study the potential for this is to check the consistency of a hypothetical observed waveform with the predicted 
higher-order terms in the general relativistic PN sequence, which depend on very few parameters (only the two 
masses, for non-spinning, quasi-circular inspirals). Another is to examine the constraints that could be placed on 
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specific alternative theories using gravitational- wave observations [2-10]. Most of these analyses have incorporated 
only the dominant effect that distinguishes the chosen theory from general relativity, such as dipole radiation or the 
wavelength-dependent propagation of a massive graviton (see, however [11]). 

Some authors have taken a different approach by proposing parametrized versions of the gravitational waveform 
model [12-15], inspired by the parametrized post-Newtonian (PPN) formalism used for solar-system experiments, 
and analysing the bounds that could be placed on those theory-dependent parameters by various gravitational-wave 
observations. Yet the authors of these frameworks were limited by the fact that for many alternative theories of 
gravity, only the leading terms in the waveform model have been derived. 

In addition, the existing parametrizations of the gravitational waveform make the implicit assumption that the 
gravitational wave signal during the inspiral depends only on the masses of the orbiting compact bodies (in the 
spinless case), and not on their internal structure. This is true in general relativity, which satisfies the Strong 
Equivalence Principle, but is known to be violated by almost every alternative theory that has ever been studied. In 
scalar-tensor theory, for example, the internal gravitational binding energy of neutron stars has a definite effect on the 
motion and gravitational- wave emission, and since the binding energy can amount to as much as 20 percent of the total 
mass-energy of the body, the effects can be significant. In order to determine the full nature of the gravitational- wave 
signal in an alternative theory of gravity, the strong internal gravity of each body must be accounted for somehow, 
even in a PN expansion. 

To make the situation even more interesting, binary black holes play a special role within the scalar-tensor class 
of alternative theories. Based on evidence from a 1972 theorem by Hawking [16], together with known results from 
first-post-Newtonian theory, it is likely that in a broad class of scalar-tensor theories, binary black hole motion and 
gravitational radiation emission are observationally indistinguishable from their GR counterparts. This conjecture will 
be discussed in more detail later in this paper. 

Scalar-tensor gravity is the most popular and well-motivated class of alternative theories to general relativity. Apart 
from the long history of such theories, dating back more than 50 years to Jordan, Fierz, Brans and Dicke [17], scalar- 
tensor gravity has been postulated as a possible low-energy limit of string theory. In addition, a wide class of so-called 
f(R) theories, designed to provide an alternative explanation for the acceleration of the universe to the conventional 
dark-energy model, can be recast into the form of a scalar-tensor theory (for reviews, see [18, 19]). 

Measurements in the solar system and in binary pulsar systems already place strong constraints on key parameters 
of such theories, notably the coupling parameter uj . Yet these tests probe only the lowest-order, first post-Newtonian 
limit of these theories, some aspects of their strong-field regimes (related to the strong internal gravity of the neutron 
stars in binary pulsars) and the lowest-order, dipolar aspects of gravitational radiation damping. 

These considerations have motivated us to develop the full equations of motion and gravitational waveform for 
compact bodies in a class of scalar-tensor theories to a high order in the PN sequence; this is the first in a projected 
series of papers aiming to treat this problem in full. 

It should be acknowledged that we do not expect any big surprises. Damour and Esposito-Farese [20] have shown 
on general grounds that the available constraints on the scalar-tensor coupling constant loq derived from solar-system 
experiments imply that scalar-tensor differences from GR will be small to essentially all PN orders, except for certain 
regions of scalar-tensor theory space where non-linear effects inside neutron stars, called "spontaneous scalarization" , 
can occur (for a recent example, see [21]). It is therefore unlikely that we will be able to point to a qualitatively new 
test of scalar-tensor gravity to be performed with gravitational waves. 

Nevertheless we expect to provide a complete and consistent waveform model to an order in the PN approximation 
comparable to the best models from GR. With this model it will be possible to carry out parameter estimation 
analyses for gravitational waves from binary inspiral, and to compare the bounds with those from earlier work that 
cither confined attention to the leading dipole term, such as [6], or assumed extreme mass ratios, such as [11]. 

We will use a version of the formalism of "post-Minkowskian" theory, which has proven to be very powerful for 
deriving the equations of motion and gravitational-wave signal to high post-Newtonian orders in GR. The specific 
version is known as Direct Integration of the Relaxed Einstein Equations (DIRE), based on a framework originally 
developed by Epstein and Wagoner [22], extended by Will, Wiseman and Pati [23-26], and applied to numerous 
problems in post-Newtonian gravity [27-32]. DIRE is a self-contained approach in which the Einstein equations are 
cast into their "relaxed" form of a flat-spacetime wave equation together with a harmonic gauge condition, and are 
solved formally as a retarded integral over the past null cone of the field point. The "inner", or near-zone part of 
this integral within a sphere of radius A, a gravitational wavelength, is approximated in a slow-motion expansion 
using standard techniques; the "outer" part, extending over the radiation zone, is evaluated using a null integration 
variable. 

DIRE is rather easily adapted to scalar-tensory theory, so that the same methods that have been worked out for GR 
can be applied here. It is possible that many other theories that generalize the standard action of general relativity 
in four spacetime dimensions by adding various fields could be cast in a similar form, permitting a systematic study 
of their predictions for compact binary inspiral beyond the lowest order in the PN approximation. Indeed another 
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motivation for this work is to lay out a template for possible extensions to other theories of gravity, such as the 
Einstein- Aether theory [33] or TeVeS [34]. 

Specifically, the theories we address here are described by the action 

S = (167T)- 1 J [<j>R - 4>- l u(<t>)9 ap d*<t>dp<t>] V^d^x + S m (m,g a p) , (1.1) 

where R is the Ricci scalar of the spacetime metric g a fi, $ is the scalar field, of which u is a function. Throughout, 
we use the so-called "metric" or "Jordan" representation, in which the matter action S m involves the matter fields m 
and the metric only; cf) does not couple directly to the matter (see [35] for example, for a representation of this class 
of theories in the so-called "Einstein" representation). We exclude the possibility of a potential or mass for the scalar 
field. 

In order to incorporate the internal gravity of compact, self-gravitating bodies, we adopt an approach pioneered 
by Eardley [36], based in part on general arguments dating back to Robert Dicke, in which one treats the matter 
energy-momentum tensor as a sum of delta functions located at the position of each body, but assumes that the mass 
of each body is a function Ma (4>) of the scalar field. This reflects the fact that the gravitational binding energy of the 
body is controlled by the value of the gravitational constant, which is directly related to the value of the background 
scalar field in which the body finds itself. Consequently, the matter action will have an effective dependence on <j>, 
and as a result the field equations will depend on the "sensitivity" of the mass of each body to variations in the scalar 
field, holding the total number of baryons fixed. The sensitivity of body A is defined by 

•*-(™)- 

For neutron stars, the sensitivity depends on the mass and equation of state of the star and is typically of order 0.2; 
in the weak- field limit, sa is proportional to the Newtonian self-gravitational energy per unit mass of the body. From 
the theorem of Hawking, for stationary black holes, it is known that sbh = 1/2. 

This paper reports the results of a calculation of the explicit equations of motion for binary systems of non-spinning 
compact bodies, through 2.5PN order, that is, to order (v/c) 5 beyond Newtonian theory. The post-Newtonian 
corrections at 1PN and 2PN orders are conservative; we obtain from them expressions for the conserved total energy 
and linear momentum, and obtain the 2-body Lagrangian from which they can be derived. There are also terms 
in the equations of motion at 1.5PN and 2.5PN orders. These are gravitational-radiation reaction terms. Terms at 
1.5PN order do not occur in general relativity, but in scalar-tensor theories with compact bodies, they are the result 
of the emission of dipole gravitational radiation. At 2.5PN order, one finds the analogue of the general relativistic 
quadrupole radiation, together with PN correction effects related to monopole and dipole radiation. 

Not surprisingly the expressions for these quantities are complicated, much more so than their counterparts in 
general relativity. On the other hand, they depend on a relatively small number of parameters, related to the value 
of ui(<p) far from the system, where (j> = (fio, along with its derivatives with respect to tp = 4>/4>o, and the sensitivities 
Si and s 2 of the two bodies, and their derivatives with respect to <f>. The parameters and their definitions are shown 
in Table I. 

At Newtonian order, the "bare" gravitational coupling constant G is related to the asymptotic value of the scalar 
field, but for two-body systems of compact objects, the coupling is given by the combination Ga, where 

_ 3 + 2^ (l-2 Sl )(l-2s 2 ) 
"~ 4 + 2^ 4 + 2c.o ' 1 j 

where cj = w(<fo)- At 1PN order there are two body-dependent parameters, 7 and (3a, A = 1,2 (see Table I for 
definitions of the parameters). For non-compact objects, where s J 4 < ^l ) 7 = 7 — 1 and Pa = P — 1 ; where 7 and ft 
are precisely the PPN parameters for scalar-tensor theory, as listed, for example in [37] . At 2PN order, there are two 
additional parameters 8a and \A- Most of the parameters in Table I can be related directly to parameters defined 
in [20, 35]. 

Here we will quote the bottom-line result: the two-body equation of motion, expressed in relative coordinates, 
x = xi — X2, through 2PN order. This equation is ready-to-use, for example in calculating time derivatives of 
radiative multipole moments in determining the gravitational-wave signal, which will be the subject of the second 
paper in this series. The equation has the form 

d 2 x Gam Gam r . „ ,-, 

-^2 = —n+— 2 —[n(Ap N + A 2 p N )+rw(Bp N +B 2PN )\ 

8 (Gam) 2 r . /A . ,_, Nl ,„ A . 

+-=r\- 3-^- [rn(Ai. 5Pj v + A 2 . 5PN ) - v(Bi. 5 pjv + B 2 . 5 pn)\ , (1-4) 
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TABLE I: Parameters used in the equations of motion 



Parameter 


Definition 


Parameter Definition 






Scalar-tensor parameters 


Equation of motion parameters 






G 


<A7 1 (4 + 2u; () )/(3 + 2a;o) 


Newtonian 






C 


1/(4 + 2uj ) 


a l-C + C(l-2si)(l-2s 2 ) 






Ai 


(do;/^)oC 2 /(l - 


post-Newtonian 






A 2 


(dWKVu - 


7 -2a _1 C(l - 2si)(l - 2s 2 ) 






Sensitivities 




/3i a" 2 C(l - 2s 2 ) 2 (Ai(l - 2si) + 2^) 






SA 


[d In M A {4>)/d In 0] o 


& a~ 2 <T(l — 2si ~) 2 CAi fl - 2sc1 + 2Cs'->) 






1 
SA 


[d 2 lnM A (0)/dln0 2 ]o 


2nd post-Newtonian 






s'k 


[d 3 lnM A (0)/dln0 3 ] o 


51 a- 2 C(l-C)(l-2si) 2 

5 2 a- 2 C(l-C)(l-2s 2 ) 2 










Xi a- 3 C(l-2s 2 ) 3 [(A 2 -4A 2 + CAi)(l- 
X2 a- 3 C(l-2 Sl ) 3 [(A 2 -4A 2 + CAi)(l- 


2si) 


-6CA 1 s' 1 + 2C 2 si'] 






2s 2 ) 


-6CA lS2 + 2C 2 s 2 '] 



where r = |x|, n = x/r, m = mi + 7712, 77 = mi"i2/™ 2 , v = vi — V2, and r = dr/dt. We use units in which c = 1. The 
leading term is Newtonian gravity. The next group of terms are the conservative terms, of integer PN order, while 
the final group are dissipative radiation-reaction terms, of half-odd-integer PN order. The coefficients A and B are 
given explicitly by 



Apn 
Bpn 
A2PN 



B 



2PN 



Al.SPN 
B\.5PN 



3 Cjoltcl 
-(1 + 3?7 + 7 )u 2 + -r]f 2 + 2(2 + 77 + 7 + j3+ - #-)— — , 
2 r 

2(2-77 + 7), 

—77(3 - 477 + 7)7/ + 1 [77(13 - 477 + 4 7 ) - 4(1 - 4t7)/3+ + 4i/>(l - 3r7)/3_] v ^^H _ 1^(1 _ 3^ 



(1.5a) 



+ -77(3-4r/ + 7)7jV + 



2 + 25/7 + 2t7 2 + 2(1 + 977)7 + - 7 2 - 477(3^+ + 28+ + 2#_ 



9 + ^77 + (9 + 877)7 + 1(9 - 2t7)7 2 + (8 + 157? + 4 7 )^+ - ^(8 + 777 + 4 7 )^_ 



+ (1 - 2r,)(S+ - 2x+) + V(^- + 2x-) - 247? 

1 77(15 + 477 + 87)7/ - |t7(3 + 2t] + 2 7 )f 2 
1 



P1P2 



// Gam ^ ' 



- [4 + 4I77 + 87/ 2 + 4(1 + 777)7 + 7 2 - 87/(2/3+ - ^/3_) + 45+ + 4#_ 



Gam 



, Gam 



(1.5b) 



(1.5c) 



The expressions for the 2.5PN coefficients are lengthy and will be displayed in a later section. Here the subscripts 
"+" and "— " on various parameters denote sums and differences, so that, for a chosen parameter wc define 



T - = l(n-r 2 ) 



(1.6) 



The quantity <S_ and its companion S+ (used later) are defined by 



S- = -aT 1/2 ( S i-.s 2 ), 
S+ = a- 1 / 2 (l-s 1 -s 2 ), 



(1.7) 
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(the significance of these definitions of S± will become clear in Sec. VII A), and tp is defined by 



^=^^=±^1^. (1.8) 

TTl\ -\- 777-2 

Several things are worth noting about these equations (and indeed about all the two-body equations shown later 
in this paper). In the general relativistic limit lu — > oo, or £ — > 0, the equations (including the 2.5PN terms) reduce 
to those of general relativity, as determined by many authors [26, 38-43]. At 1PN order, the equations agree with the 
standard scalar-tensor equations, both for weakly self-gravitating bodies in the general class of theories [44] (shown 
within the PPN framework in Sec. 6.2 and 7.3 of [37]), and for arbitrarily compact bodies in pure Brans-Dicke theory 
(as displayed in Sec. 11.2 of [37]). 

Although a number of authors have obtained partial results in scalar-tensor theory at 2PN order, notably the 
metric sufficient to study light deflection at 2PN order [45, 46], and the generic structure of the 2PN Lagrangian for 
N compact bodies [20], our explicit formulae for the 2PN and 2.5PN contributions to the two-compact-body equations 
of motion are new. 

The energy loss that results from the 1.5 PN and 2.5 PN terms in the equations of motion is in complete agreement 
with the energy flux calculated to the corresponding order by Damour and Esposito-Farese [35] . 

The other interesting limit is that in which both bodies are black holes. Assuming that Hawking's result that 
sbh = 1/2 applies equally for binary black holes as for isolated black holes, we find that the parameters 7, /3 a, 5 a and 
Xa all vanish, and a = 1 — ( = (3 + 2cj )/(4 + 2aj ). In this case the equations reduce identically to those of general 
relativity through 2.5PN order, with Gam a replacing of Guia for each body. In other words, if each mass is rescaled 
by (4 + 2wo)/(3 + 2wo), the scalar-tensor equations of motion for binary black holes, including the 2.5PN terms, 
become identical to those in general relativity. Again this applies to all the equations of motion and related quantites 
(total energy, Lagrangian), whether for the individual bodies or for the relative motion. Since the masses of bodies 
in binary systems are measured purely via the Keplerian dynamics of the system, the rescaling is unmeasurable, and 
therefore, the dynamics of binary black holes in this class of theories is observationally indistinguishable from the 
dynamics in general relativity. Assuming, as we believe will be the case, that this is also true for the gravitational 
wave emission, the conclusion is that gravitational-wave observations of binary black hole systems will be unable to 
distinguish between these two theories. 

If only one member of the binary system is a black hole, then a = 1 — £, and 7 — (3a = 0, so that even at 1PN order, 
the equations of motion are identical to those of general relativity, after rescaling each mass. Only at 1.5PN order 
and above do differences between the two theories occur for the mixed binary system, because of the non-vanishing 
of iS_ in the dipole radiation reaction term, and the non-vanishing of <5i (if body 1 is the neutron star) in the 2PN 
terms. However, in this case all the deviations from general relativity depend on a single parameter Q, given by 

Q^((l-()- 1 (l-2s 1 f, (1.9) 

where si is the sensitivity of the neutron star. In particular, all reference to the parameters Ai and A2 disappears, 
and the motion through 2.5PN order is identical to that predicted by pure Brans-Dicke theory. If this conclusion 
holds true for the gravitational-wave emission, then gravitational-wave obervations of mixed black-hole neutron-star 
binaries will be unable to distinguish between Brans-Dicke theory and its generalizations. The only caveat is that, for 
a given neutron star, generalized scalar-tensor theories can predict very different values of its un-rescaled mass and 
its sensitivity from those predicted by pure Brans-Dicke. 

The remainder of this paper provides details. In Sec. II, we derive the "relaxed field equations" and the associated 
formal equations of motion in scalar-tensor theories, and write down the formal solutions for the gravitational and 
scalar fields in terms of solutions of the flat spacetime wave equation. In Sec. Ill we describe the formal structure 
of the fields in the near zone, and in Sec. IV we obtain formal solutions for the fields through 2.5PN order in terms 
of Poisson-like potentials and time derivatives of system multipole moments. In Sec. V we introduce the Eardley 
approach for characterizing the compact bodies, rewrite all equations in terms of a "conserved" density in which the 
masses of each body are constant, and arrive at the equations of motion expressed in terms of the redefined potentials. 
In Sec. VI we apply the methods of [26] to express the equations of motion explicitly in terms of masses, positions 
and velocities of each compact body in a two-body system. We obtain the 2-body Lagrangian, the conserved total 
energy and linear momentum, the relative effective one-body equations of motion, and the rate of energy loss due to 
gravitational-radiation reaction. Section VII presents a detailed discussion of the results. 
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II. THE RELAXED FIELD EQUATIONS IN SCALAR-TENSOR THEORY 
A. Field equations and equations of motion 

We begin by recasting the field equations of scalar-tensor theory into a form that parallels as closely as possible 
the "relaxed Einstein equations" used to develop post-Minkowskian and post-Newtonian theory in general relativity. 
The original field equations of scalar-tensor theory as derived from the action of Eq. (1.1) take the form 

Gp.v = ^-T„ v + ^ (4> tli <t>, v ~ \g^<t>,x^ + ^ {4>;n» - 9n^g<t>) , (2.1a) 

where T^ v is the stress-energy tensor of matter and non-gravitational fields, G M „ is the Einstein tensor constructed 
from the physical metric (/ M „, <f> is the scalar field, <jj(4>) is a coupling function, d fl denotes the scalar d'Alcmbcrtian 
with respect to the metric, and commas and semicolons denote ordinary and covariant derivatives, respectively. We 
work throughout in the metric or "Jordan" representation of the theory, in contrast to the "Einstein" representation 
used, for example in [35]. 

Normally, such as for a perfect-fluid source, the matter stress-energy tensor depends only on the matter field 
variables and the physical metric <? M „, not on the scalar field, and accordingly the term dT/d<j> does not appear in 
the field equations. But in dealing with a system of self-gravitating bodies, we will adopt an approach pioneered by 
Eardley [36] . Because <p controls the local value of the gravitational constant in and near each body in this class of 
theories, the total mass of each body, including its self-gravitational binding energy, may depend on the scalar field. 
Thus, as long as each body can be regarded as being in stationary equilibrium during its motion, Eardley proposed 
letting each mass be a function of <fi, namely Ma{4>)- With this assumption, T^ v takes the form 

T^(x a ) = {-g)- 1/2 Y. I ' drMA^uW^xlir) - x a ) 

A J 

= (-g)- 1/2 Y, M ^ u >A(u A r 1 S 3 ^-^A), (2.2) 

A 

where r is proper time measured along the world line of body A and u A is its four- velocity. The indirect coupling of 
4> to matter via the binding energy is responsible for the term dT/dcj) in the field equations. 
From the Bianchi identity applied to Eq. (2.1a), the equation of motion is 

T^;v = , (2.3) 

with the right-hand-side vanishing in the perfect-fluid case. From the compact body form of T^ v in Eq. (2.2), it can 
then be shown that the equation of motion for each compact body takes the modified geodesic form 

u v V„{M A {<j>)u») = -^&» , (2.4) 
acp 

or in terms of coordinate time and ordinary velocities v a , 

£ + r>v» - r>"„V . -j^^-' - *V> . (2-5) 

These equations of motion could also be derived directly from the effective matter action, S m — J2a Ia^ a ^^ t - 
Equation (2.2) can equally well be taken to describe a pressureless perfect fluid (dust), simply by letting the mass of 
each particle be a constant, independent of <j>. 



B. Relaxed field equations 

To recast Eq. (2.1a) into the form of a "relaxed" Einstein equation, we make use of the following well-known 
property: defining the quantities 

jT = >/=S<r, (2.6a) 
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Rmv ji = s w> ga p_ g av g pp t (2.6b) 
it can be shown that the following is an identity, valid for any spacetime, 

a0 = (-g)(2G^ + l6nt" L v L ) , (2.7) 

where t^ L is the Landau-Lifshitz pscudotensor [see Eq. (20.22) of [47] for an explicit formula]. 

To incorporate scalar-tensor theory into this framework, we assume that, far from any isolated source, the metric 

takes its Minkowski form r/^, and that the scalar field <fi tends to a constant value <f>Q. We define the rescaled scalar 
field ip = 4>/(f>o- We next define the conformally transformed metric g^ u by 

9nv = ipg^ , (2.8) 

and the gravitational field by the equation 

= yf^gcT = rT - h"" . (2.9) 
From Eq. (2.8) it can be shown that this is equivalent to 

i^^ip- 1 ^ -hT). (2.10) 

We now impose the "Lorentz" gauge condition 

W,„ = 0, (2.11) 

which is equivalent to 

= -<p- 2 <p, v (*r - h * u ) ■ ( 2 - 12 ) 

Substituting Eqs. (2.1a), (2.1b), (2.9) and (2.11) into (2.7), we can recast the field equation (2.1a) into the form 

n v h» v = -Wttt^ , (2.13) 
where D v is the flat spacetime d'Alembertian with respect to rj^ Ul and where 

16ttt^ = 16tt(-. 9 )^T^ + A"" + Ag" , (2.14) 

where 

ST = 16n[(-g)t^ L }(3n + h^, h^, a -h a0 h^, a0 , (2.15a) 

AgT ee £^l<p M - , (2.15b) 

where the notation [(— g)tj* / L \(g' J ' 1 ') denotes that the Landau-Lifshitz piece should be calculated using only q, in other 

words, exactly as in general relativity, except using the conformal metric, rather than the physical metric. The scalar 
field equation can also be rewritten in terms of a flat-spacetime wave equation, of the form 

□„¥> = -8itt s , (2.16) 

where 



3 + 2uj v a (j) V dip J 8tt 



J d_ 

167r dip 



hi 



2u 



f 2 



VMQ a0 ■ (2-17) 



In principle, Eqs. (2.6a) and (2.9) can be combined to give g^ v in terms of ip and W v , although in practice, we will 
express it as a PN expansion. The final result will be the relaxed field equations (2.13) - (2.17) expressed entirely 
in terms of h^, ip, and the matter variables. The next task will be to solve these equations iteratively in a post- 
Newtonian expansion in the near-zone. Formally the solutions of these wave equations can be expressed using the 
standard retarded Green function, in the form 

J |x-x'| 

¥>(*,*) = 2/^^j^V, (2.18) 

where the integration is over the past flat spacetime null cone of the field point (t, x). We will expand these integrals 
in the near-zone, and incorporate a slow-motion, weak- field expansion in terms of a small parameter e <~ v 2 <~ m/r; 
the strong-field internal gravity effects will be encoded in the functions Ma (</>)■ 
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III. FORMAL STRUCTURE OF THE NEAR-ZONE FIELDS 

We follow [25] (hereafter referred to as PWI) by defining a simplified notation for the field W and the scalar field 



N 


EE h m 




K l 


EE h 0i 






EE h 13 


~ 0{e 2 ) , 


B 


EE h U ■ 


i 




= V- 


l~0(e), 



o( £ 2 ), 

(3-1) 

where we show the leading order dependence on e in the near zone. To obtain the equations of motion to 2.5PN 
order, we need to determine the components of the physical metric and <p to the following orders: goo to 0(e 7 / 2 ), goi 
to 0(e 3 ) , Qij to 0(e 5 / 2 ), and ip to 0(e 7 / 2 ). From the definitions (2.6a) and (2.9), one can invert to find g^ v in terms 
of W and ip to the appropriate order in e, as in PWI, Eq. (4.2). Expanding to the required order, we find, 

.9oo = -1 + QiV + *) + Qs - ^iV 2 - iiVtf - * 2 

+ ( —N 3 - -NB + -K j R3 + -NH - -B^> + -N^ 2 + * 3 
V 16 4 2 8 2 2 

+0(e 4 ) , (3.2a) 
5w = -K* + QiV + K l + 0{e 7 ' 2 ) , (3.2b) 

9ij = S ij |l + QiV - - QiV 2 + Is + I TV* - ^ | + S« 

+0(e 3 ) , (3.2c) 
(-g) = 1 + (N - 4*) - (B + 4N<$ - 10* 2 ) + C>(e 3 ) . (3. 2d) 

In Eqs. (3.2) we do not distinguish between covariant and contravariant components of quantities such as K % or B n ' , 
since their indices are assumed to be raised or lowered using the Minkowski metric, whose spatial components are Sij. 

We now define a set of provisional "densities" following the convention of Blanchet and Damour [48] , but adding a 
separate density for the scalar field equation: 

a ee T 00 +T U , 
a' ee T°\ 

a s ee -T + 2pdT/dip. (3.3) 

The second contribution to a s will be non-zero only in the case where our system consists of gravitationally bound 
bodies, whose internal structure could depend on the environmental value of <p. 

Because of the way we have formulated the relaxed scalar-tensor equations, the quantity A M!/ has exactly the same 
form as in PWI, Eq. (4.4). To the 2PN order needed for our work, we have 

A 00 = --(ViV) 2 + \-N 2 -NN- 2N' k K k + + K 1 - 1 ) 

8 8 2 

+k i N> i - B ij N' ij + J ViV • VB + ^N(\7N) 2 J + 0(pe 3 ) , (3.4a) 

A 0i = I^N' k (K k > 1 - K l ' k ) + ^NN'^ + 0(pe 5 / 2 ) , (3.4b) 

A ij = h^N'iN'i - i^'(ViV) 2 } + 1 2 A h ■ ' A ' ■* - A ' ' A' ' - K i ' k K j ' k + 2N^K j '> + ^N'^B'^ 

--N(N' l N J - -S ij (VN) 2 ) - 5 ij (K l ' k K^ + N< k K k + -N 2 + -VN ■ VB)\ + 0{pe 3 ) , (3.4c) 
2 2 8 4 J 



9 



A" = -^(ViV) 2 + ^K l ' k K^ - N' k K k - ^\7N • V5 - ^N 2 + ^N(\7N) 2 ^ + 0(pe 3 ) . (3.4d) 



To the required order, the scalar stress-energy pseudotensor is given by 



A oo = 3 + 2wo (w)2 + 3 + 2^, ( N(y ^2 _ 2 (i- -4-") *(V*) 2 + vj/ 2 } 

+0(pe 3 ) , (3.5a) 

- -(3 + 2w )**^ + O(pe 5 / 2 ), (3.5b) 

A| = (3 + 2w ) jV 1 ^ - i^'(V*) 2 



-(3 + 2w ) 2 1- 



3 + 2w / 



+0(pe 3 ), (3.5c) 



AS = -^^(V*) 2 + (3 + 2c ) { (l - ^|^) *(W) 2 + + 0(p, 3 ) , (3.5d) 



where co' = (dw/dip)o- 

The near- zone expansions of the fields N, K l , B l i and are then given by 



N H = 4e f T ™ {t ^'} d?x' + 2e 2 d 2 t [ r 00 (t, x')|x - x'|dV - \e>' 2 l kk (t) + l -e 3 dt [ T 00 (t, x')|x - x'| W 
l x ~ x I Jm 6 J M 

r (5) (5) (5) i 

-^e 7/2 ^ (4^ + 2r 2 8 kl ) X k \t) -Ax k X ku {t) +2 kkll (t) \ + N aM + 0(e 4 ) , (3.6a) 

XV - 4e 3 / 2 y \ x *,\ d3x ' + 2e5/2 % J T 0l (t, X ')\x- X '\d 3 x' + ^e 3 hx k F k (t)~F kk (t)+2e mlk J mk (t)\ 

+Kl M + 0(6 7 / 2 ) , (3.6b) 

B % = 4e 2 / T [ 3 ( *' ^ riV - 2e 5 / 2 (t) +2e 3 d 2 / r y (t,x')|x - x'|rfV 
./.m l x ~ x I Jm 

1 f (5) (5) (4) (3) 1 

__ e 7/2 J 3r 2 -pj^ _ 2a ,fe pJ k^ _g x k e mk(i jm\j)^ + g Jtfijkk^) L + B ^ + { e i ) , (3.6c) 

*N = 2e / Z*£i^dV - 2e 3 / 2 M s + e 2 d 2 / r a (t, x')|x - x'|dV 

l x — x I 

1 / (3) (3) (3) \ 1 /■ 

__ e 5/2 / r 2 Ma(t) _ 2a J22(t) +Xf (t)J +-e 3 9 4 y^r s (^x')|x-x'|W 

i f (5) (5) (5) (5) (5) ^ 

-— e 7/2 ^ r 4 M a (t) -4rV V s {t) +{Ax kl + 2r 2 <^) 2 kl (t) -4x k lf l {t) + 2* fc "(t) \ + 0(e 4 ) , (3.6d) 
where we define the moments of the system by 



lQ EE 


/ T ° X Qd 3 X, 


(3.7a) 




JM 




jiQ = 


e lab f T 0b x aQ d 3 x, 


(3.7b) 




JM 




AT iQ EE 


[ T^X Q d 3 X , 


(3.7c) 




JM 




2? EE 


/ T s x Q d 3 x , 


(3.7d) 








M s EE 


/ r s d 3 x . 


(3.7e) 




JM 
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The index Q is a multi-index, such that x® denotes x 11 . . . x lq . The integrals are taken over a constant time 
hypersurface M. at time t out to a radius 1Z, which represents the boundary between the near zone and the far zone. 
The structure of the expansions for N/s, Kfc and B 1 ^ is identical to the structure in PWI because the source r Miy 
satisfies the conservation law ; „ = 0, a consequence of the Lorentz gauge condition. However, no such explicit 
conservation law applies to r s ; nevertheless, in a post- Newtonian expansion, we will be able to show, for example, that 
the term e 3 / 2 M s actually vanishes to lowest PN order, and thus contributes only beginning at e 5 / 2 order; the other 
terms involving time derivatives of M s will also be boosted to one higher PN order. The time derivatives of the dipolc 
moments T{ do not vanish in general; this is related to the well-known phenomenon of dipole gravitational radiation 
that can occur in scalar-tensor theories. The boundary terms Nqj^, Kq M and B % Jj M can be found in Appendix C of 
PWI, but they will play no role in our analysis. As in PWI, we will discard all terms that depend on the radius 1Z of 
the near-zone; these necessarily cancel against terms that arise from integrating over the remainder of the past null 
cone; those "outer" integrals can be shown to make no contribution to the near zone metric to the PN order at which 
we are working. 

In the near zone, the potentials are Poisson-like potentials and their generalizations. Most were defined in [25], but 
we will need to define additional potentials associated with the scalar field. For a source /, we define the Poisson 
potential to be 



P ^ S tJ m W^- 3 "'' V*P(/) = -/. (3.8) 
We also define potentials based on the "densities" a, a 1 and and a s constructed from T a/3 and from T — 2LpdT/dip, 

] -d 3 x' = P(47rcr/) , (3.9a) 



n.n = 


J 

J M 




j 




s 




j 



a(t,x>)f(t,x') j3 , 



«f« = «, (3.90 
"■ ( '' X,)/ !^<iV = P(4™,/), (3.9d) 



X — X' 



along with the superpotentials 

X(f) = f a(t,x')/(t,x')|x-x'|dV, (3.10a) 

JM 

Y(f) = f .(t,x')/(t,x')|x-x'|¥x', (3.10b) 

JM 

and their obvious counterparts X 1 , X Sl and so on. A number of potentials occur sufficiently frequently in the PN 
expansion that is it useful to define them specifically. There are the the "Newtonian" potentials, 

JM l x — x I 

U s = [ ^l^W-P^H^a). (3-llb) 

JM l x — x I 

The potentials needed for the post-Newtonian limit are: 

V 1 = S l (l), $[ j = $x = S"(l), <f>l = Y, s {v 2 ), 

$ 2 = £(co, a>5 = s s (c/) 7 $ 2s = E(u a ) , n s = ^s(u s ), 

X = X(l), X s = X s (l). (3.12) 



Useful 2PN potentials include: 



$| = S(^), Y = Y(1), 

X^X^l), X 1 = X ii (l), 
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X 2 = X(U) , 
P l 2 3 = PiU'V'i) , 



X 2s 
X 2s 



X(U S ), 
X S (U S ), 



P 2 EE P 2 " = $ 2 - -U 2 , 



pij - 




P 2s 


= PS, = *L - 2^ 


Gl 




Gis 




G 2 


= p(uti), 


G 2s 




G 3 = 


-P{U- k V k ), 


G3 S 


= -P(U*V k ) , 


G4 = 




G 5 ; 


= -P(l> fc C/' fc ) , 


G 6 = 


/': (•''<!>,•';. 


G6s 


EE P(Uf^) , 



G 4 7 ee P(U' k V k ' 1 ) + \P{U^U) , 

H EE Pip™ P%) , P s 

//~ = p(u™p?) , //; 



p^p^j, 
p(^p^). 



(3.13) 



IV. EXPANSION OF NEAR-ZONE FIELDS TO 2.5PN ORDER 



In evaluating the contributions at each order, we shall use the following notation, 

N = eiNo + eN 1 + e 3 / 2 N 1 . 5 + e 2 N2 + e 5 / 2 N 2 . 5 ) + 0(e 4 ), (4.1a) 

K i = e 3/2 iKi+eKi + e 3/2 Ki 5 ) + Oie 7/2 ), (4.1b) 

B = e 2 (B l + e^B^ + eB 2 + e 3 / 2 P 2 . 5 ) + 0(e 4 ) , (4.1c) 

B*j = e 2 (P^'+e 1 / 2 P^ 5 ) + 0(e 3 ), (4.1d) 

* = e(*o + e 1/2 *o.5 + e*i + e 3/2 *i.5 + e 2 *2 + e 5/2 *2.5)+G(e 4 ), (4.1e) 



where the subscript on each term indicates the level (1PN, 2PN, 2.5PN, etc.) of its leading contribution to the 
equations of motion. 



A. Newtonian, 1PN and 1.5PN solutions 

At lowest order in the PN expansion, we only need to evaluate r on = (— g)T no + O(pe) = a + O(pe) (recall that 
a %l <~ ecr), and t s = a s /[(j) (3 + 2wo)], where uj = lo{4>q). Since both densities have compact support, the outer 
integrals vanish, and we find 

A^o - ^, (4.2) 

<Po 

*o = ; ■ (4.3) 

0o (3 + 2w ) 

Consider the case where we are dealing with pure perfect fluids, with no compact bodies having sensitivity factors 
sa- Then to Newtonian order, a = a s , U = U s , and the metric to Newtonian order is given by the leading term in 
Eq. (3.2a), 

5oo = -1+ + ( 4 - 4 ) 

= -i + 2 4 + 2^ „ (45) 
0o(3 + 2w o ) V ; 

We therefore identify the coefficient of U in g 00 as the effective Newtonian gravitational coupling constant, G, given 
by 

G ee ^7—77-^ • (4.6) 
0o 3 + 2wo 
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However, we will not set G — 1 as is conventional in general relativity, in order to highlight the fact that it is an 
effective gravitational constant linked to the asymptotic value of 0, which could, for example, vary with time as the 
universe evolves. For future use, we also define the parameters 



c 

Ai 



1 



4 + 2cj ' 
(duj/d(p)o( 
3 + 2w ' 



A consequence of these definitions is that 



= GC. (4.8) 



0o (3 + 2lo ) 

It is worth pointing out that ojq enters at Newtonian order, via the modified coupling constant G of Eq. (4.6). It is 
then clear, by virtue of the expansion u){<j>) — ijj + {du/dip)^ + (d 2 ui/dip 2 )o^ 2 /2 + . . . , that the parameter Ai will 
first contribute at 1PN order, A2 will first contribute at 2PN order, and so on. 

To this order, (-g) = 1 + 4GU(1 - Q- 8GU S ( + 0(e 2 ). Then, through PN order, the required forms for and 
t s are given by 

r 00 = G(l-C)^-« tl + G(l-C)(^U-^(VU) 2 )-G((6*U s -^(VU s ) 2 )}+0( P e 2 ), (4.9a) 

t w = G(l - CV + 0(pe 3/2 ) , (4.9b) 

r H = G(l C){^ - ^G(l - C)(VC/) 2 - i-GC(VC/ s ) 2 | + 0( P e 2 ) , (4.9c) 

t 13 = O(pe), (4.9d) 

1 

2^ 



t s = Gc|<t s + 2G(1 - ()a s U - 2G(2X 1 + ()a s U s + ^G(Ai - C)(Vt/ s ) 2 ) + O^e 2 ) . (4.!Jo) 



Substituting into Eqs. (3.6), and calculating terms through 1.5PN order (e.g. 0(e 5 / 2 ) in N), we obtain 
iVi = G(l-C)|7G(l-C)C/ 2 -4$i + 26(1-0*2+21 

-G(U 2 - 24GC$ 2s + 2GC$L} , (4-lOa) 
K\ = 4G(1-C)^ 1 , (4.10b) 
B 1 = G(1-C)|g(1-C)C/ 2 +4$i -2G(1-C)*2 

+GCf/ 2 -2GC$L}, (4-lOc) 
G(( -2G(Ax - QU 2 + 4G(1 - CW2 - 4G(Ai + 2C)$* S + X s | , (4.10d) 

(4.10e) 

(4.10f) 
(4.10g) 



*1 


= GC|-2G(A 




9 (3) 


JV1.6 


= -§ , 




(3) 


Si. 5 


= -2 I fefe (t) , 


*1.5 


= -2M,(t) + - 



9 (3) 1 (3) 

_^If(t)--lf( 



In Eq. (4.10g), we have used the fact (to be verified later) that, because of the conservation of baryon number, and 
assuming that our compact bodies have stationary internal structure, M s (t) is constant to the lowest PN order. Thus, 



13 



rather than contributing to ^0.5 as shown in Eq. (3.6d), the term —2M S contributes to similarly the term in 

\I>1.5 involving three time derivatives of M s actually contributes to *2.5- 
The physical metric to 1.5PN order is then given by 

.goo = -l + 2G(l-C)t/ + 2GC(7 s -2G 2 (l-C) 2 t/ 2 -2G 2 C(C + Ai);y s 2 
-4G 2 C(1 - ()U U s + 4G 2 C(1 - C)$2 - 12G 2 C(1 - C)$2 S 
-4G 2 C(2C + Ai)$5 a + G(l - ()X + G(X S 



9m 

9ij 



A (3) o (3) 1 (3) 

-- l kk (t) -2M s (t) + -x> V s (t) -- l kk {t) +0(e 3 ) . 

-4G(1 - QV i + G(e 5/2 ) , 

% [1 + 2G(1 - ()U - 2G(U S ] + 0(e 2 ) . 



(4.11a) 

(4.11b) 
(4.11c) 



2PN and 2.5PN solutions 



At 2PN and 2.5PN order, we obtain, from Eqs. (2.14), (3.4) and (3.5), 

= G(l-()a ij + ^G 2 {l-Q 2 (U' l U> 3 -^5 ij {VU) 2 ) 
+ 4^G 2 C(1 - Q{U?U? - \&{VU a f) + 0(pe 2 ) , 



JOi 



G(l - (V + G 2 (l - o 2 (^u + lu^vw + ^uu*) 

-G 2 C(1 - 0(^U S + ^U a U?) + 0(pe 5 / 2 ) . 



(4.12a) 



(4.12b) 



Outer integrals and boundary terms contribute nothing, so we obtain 

B% = 4G(1 - C)$l j + G 2 (l - C) 2 [4P 2 ij - ^(2$ 2 - U 2 )] + G 2 C(1 - C) [*P& - ^'(2*5. - U 2 )] , (4.13a) 
K\ = G 2 (l - 2 (8Vj - 8$ 2 + 8UV 1 + 16G 4 7 ) + 2G(1 - QX* - G 2 ((l - 0(24*3. + 4P{U,U'jj) , (4.13b) 

(3) 

-2 T j (t) , 



n 2.5 



(4) 
„k n-ik 1 



(4) 



(3) 

K l 2 5 = ^x k l lk (t) ~ T kk {t) +l e mik J mk (t) 



(4.13c) 
(4.13d) 



All solutions obtained so far must be substituted into Eqs. (2.14), (2.17), (3.4) and (3.5) to obtain r 00 , t" and r s to 
the required order, 



_oo 
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G{l-Q\a- u %% + G(l - Q{AgU - ^(V£/) 2 ) - G((6aU s - ^(V*7.) 



+G 2 (l-(f{„ 
"5 • 



7G(1 - C)^ 2 - 8$i + 2G(1 - C)*2 + 2X 



- Aa ll U 



1 

47T 



U 2 - 4UU - 8U' k V k + 2V l > J (3V J ' 1 + V 1 ' 3 ) + AV J U-' 3 - W + 8WU ■ V$i 



■ VI - G(l - C) (4VJ7 • V$ 2 + 10U(\7U) 2 + AU' t] P^ J 



+G 2 C(W) \<r 



G(6Ai - 1 + 19()U 2 - G(l - C)(24UU S + 24$ 2s + 12$* 



+2G(6A! + 1 + 110*2. - 3 ^ 



+ 6a"U s 



G(l - C) (2U{\7U S ) 2 + 4U S VU • VU S + 42VC/ • V$ 2 , + 2VU S • V$ s 2 - 4VC/ • V$5 8 - 1/' </';;) 



1 

+ -T 

+ hj 2 - 2G(Ai + 2QVU S ■ V*§ 8 + ^VU S ■ VX S 
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+G(1-C)<!<7 
+0( P e 3 ) , 



(3) 
rkk 



(3) 



(3) 
kki 



- T kk (t) +6M s (t) - 2x> xi{t) +2j*(t) 



1 (3) -, (3) 

-U'^(t)+—u^At) 



r" (7i! - ,)^"-i-G(l-C)(VC/) 2 -^GC(Vf/ s ) 2 



Ua u U-— -U 2 +4V i ' j V [i ' i] +4V j U' j + -VU -VX \ 
I 47r [ 2 2 J J 



+G 2 (l-C)^4cr M C/ 

G 2 C(l-C){6a^ s -i- 



^ 2 - G(l - ()(2VU S • - 6VU • V$ 2s ) 



1 (3) 

-^(i-0^(*) +0 ^ e3 )< 

Gc|ct s + 2G(1 - ()<r s U - 2G(2X 1 + ()a s U s + ^G(Ai - C)(VC/ S ) 2 



-0 



2(1 - C)C/ 2 - 4(2Ai + 0(^ + *5) - 12C4>2 S 



(1 - C)(4*i -X) 



+G(20A 2 - 4A 2 + 6CAi + 2(?)U 2 S + 4G(2X 1 + Q(Xi + 2C)*L - (2Ai + C)*« 



-— G 2 C<{ (1 - C)(8t/t> s + Hi 17';' + ><!•', T;") + 4(Ai - Q(U 2 - VU S ■ VX S ) 



-G(l-C) 



16(Ai - C)W S • V$J - 8(1 - CWP? - SCUf'Pg 



(4.14a) 



(4.14b) 



+16G(Ai + 2C)(Ai - C)VC/ S • V$* s - 8G(A 2 - 4A 2 + 4(A : - ( 2 )U S (VU S 



(3) 
rkk 



(3) 



(3) 



-C X fefe W +^(2Ai + C) 6M.(t) - 2x^ 22 (t) + 2j*(t) 



1 (3) -I (3) 

+ -Ct/^ I«(t) +-(Ai - QU? n{t) \ + 0{pe*) . 



(4.14c) 



Substituting into Eqs. (3.6a), (3.6c) and (3.6d) and evaluating terms through 0(e 7 / 2 ), and verifying that the outer 
integrals and surface terms make no 7\L-indcpcndent contributions, we obtain, 



1 (4) 

N 2 = G{l-Q{ - Y -2Xr + G(l-0 



-4Gi - 16G 2 + 32G 3 + 24G 4 - 16G 5 - 16G 6 



7UX - 16f/*i - W l V l - 16S($i) + S(X) + 8E<(V 

20 



+ G 2 (1-C) Z 



8C/$ 2 + — U 3 - 1QH 



+G 2 C(1 - ()< 24E"(t/ s ) - (7 S 1 S - 12E(X 5 ) + E S (X S ) - 12l 2s + X| 5 + 4G is 



+G(1-C) 



8I7$2 a - 4f7 s *l - 84C/$ 2s - 4C/C/ 2 - 12E($ 2s ) - 48£($ 2 ; ) 



+4E S ($ 2 ; ) + 4S S (C/C/ S ) - 12E(J7J7 S ) - 16ff s 



24G(Ai + 3C)E(t/ 2 ) 



+4G(Ax + 2C) 



12S($y - £,(*«,) - £ s ([/ 2 ) + C/ s d>L 



B 2 = G(l-0 2X 1 + G (l-0 



t/x + 4y 4 y l - s(x) - 8E*(^ 1 ) + i6S"(c/) 



-X 2 - 20Gi + 8G 4 + 16G 5 



} + G 2 C(i-o{ 



(4.15a) 



+ G 2 C(1 - U S X S - 24£«([/ s ) - £ s (l s ) - X s 2s + 4G ls 



15 



B2.5 

*2.5 



+G(1-C) 



4U S $ S 2 - 12C/$ 2s + 12£($ 2s ) - 4E a ($5) - 4E a (IW a ) + 12E(UU a 



+4G(Ax + 2C) 



£ s ($* s ) + £ s ([/ 2 )-;7 s <l>L 

2^ S {X) - 8E 8 (*i) - 8G 2s + 16G 3s - 8G 6s + 2X 2 S 
2G(Ax + 2C)(S S (1 S ) + X° s ) - 2G(Ai - C)^* s - 8G 2 (1 - 0(Ai + 2C) 



(4.15b) 



f 1 ( 4 ) 

Gc|-n+G(i-o 



£ a (*5) + £ a (EW a 



(A 1 -C)f/ s $L + (Ai+2C)S s (^ 



8G 2 (1-C)(A!-C)^^ 



+8G 2 (Ax +2C) 

+G 2 (1 - 2 (4S S ((7 2 ) - 8tf s ) - G 2 C(1 - 0(24S s ($ 2s ) + 8H S S ) 
-^G 2 (A 2 - 4A 2 + 4CAx - C 2 )^ 3 - 4G 2 (A 2 - 4A 2 - 5CA X - 4C 2 )S S (C/ 2 ) 



(4.15c) 



(5) 



(5) 



^2.5 - 



-(2x kl + r 2 S kl ) l kl {t) + —x k X ku (t) 

15 15 

(3) 



(5) 



rkkU 



30 



(3) 



(*) +G(1 - 



(3) 



16, 



(3) 
rkk 



(3) 
fcZ n~kl 1 



U T kk {t) -AX- kl T kl {t) 



+2AUM s (t) - 8(x k U - X> k ) l k (t) +4U l kk (t) —£ x i k X s(*) 



(5) 
2 -rkk. 



(5) 

k s-rkll I 



(4) 



(3) 



(3) 
,k <rk 1 



~r 2 l kk (t) + ^x k l ku (t) + -x k t mk ° J mj (t) ~ M kku (t) +^G(1 - ()Xf l k (t) , 



(4.15d) 
(4.15c) 



I (5) 1 (5) , (5) (5) .. (3) 

■-(2x kl + r 2 5 kl ) l kl (t) +-x k l kl \t) -- l kkll (t) + -r 2 x k l k (i) --r 2 M s (t) 



(5) 



(5) 



(5) 



+GC 



(3) 



(3) 



jU, l kk (t) -2Xfl kl {t) 



4 ( 3 ) 

-G(A 1+ 20^I s fe (t) 



+ -G{2X l + QU s 



(5) (3) ■ 

6M s (t)~2x k l k (t) +l kk (i) 



(4.15f) 



V. PN EXPANSION OF THE MATTER SOURCE 



A. Energy momentum tensor and the conserved density 



We now must expand the effective energy-momentum tensor, Eq. (2.2) in a PN expansion to the required order, 
including the <j> dependence of the masses Ma- We first expand Ma{4>) about the asymptotic value (f>o- 



M A {<j ) )=M A0 + 8<j>[^^j 
We then define the dimensionless "sensitivities" 

SA = 



s'a 



1 „ l0 fd 2 M A 



dlnM A (<j>) \ 
d]n<f> J ' 

d 2 \nM A {4>) \ 
d(ln«£) 2 ) - 

d 3 In M A {4>Y 
d(\ncj)) 3 



1 3 / W1V1A 



6^ ^ 



+ ... . 



(5.1) 



(5.2) 



Note that the definition of s' A used in [37] and [49] has the opposite sign from our definition. Recalling that 
4>o (1 + ^)we can write 



M A {4>) = m A 



I + sa^ + ^a + s'a-sa)^ 2 
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+ g(4 + 3sW - Zs'a + 4 - 34 + 2^)* 3 + 0(* 4 ) 

= m A [l + S(s A ;^)} , (5.3) 

where we define the constant mass for each body rriA = Mao- 

In general relativity, neglecting pressure, the stress energy tensor can be written as (PW II, Eq. (2.12)) 

= p*{-g)- x l' 1 u> 1 u v /u° , (5.4) 

where p* is identified as the "baryonic", or "conserved" mass density, p* = mn^/^g u° , where n is the number density 
of baryons, and m is the rest mass per baryon. It satisfies an exact continuity equation dp* /dt + V • (p*v) = 0, and 
implies that the baryonic mass of any isolated body is constant. Here we identify the "baryons" as our compact point 
masses with constant mass tua, so that 



p* = ^2m A 5 3 (x- x A ) , 



(5.5) 



Thus, we can rewrite Eq. (2.2) in the form 

= p'f-f 1 / 2 /// [1 + S(s; *)] , 



(5.6) 



where p* is given by Eq. (5.5), and where we have substituted u M = w v M , with = dx^/dt — (l,v) being the 
ordinary velocity. We have dropped the subscript from the variable s in S because it will be assigned a label A 
wherever the delta function that is implicit in p* corresponds to body A. Thus, we arrive at a conversion from the 
cr-densities of Eq. (3.3) to p* , given by 



(7 = 

a 1 = 
a 11 = 



p^-gr^il + v^ll+S^n , 

p^-gy^Wii + s^n , 

p^-g)- 1 ' 2 ^^^ [!+$(«;*)] . 



To convert er s , recall that 



T = 



9^ v 

and that ip = 1 + d/d<p — d/d^. Consequently 



-T + 2<p— 

dip 



= p'i-g)- 1 '^ )- 1 
= p*(-g)- 1 / 2 (u a )- 1 [(l-2 S )+S s (s^)} 



BS 

1 + S - 2 ( 1 + ^d* 



Defining 



a s = s 2 + s' - -s , 

a s ' = s" + 2ss' - is' , 
b s = aj - a s + sa s , 



(5.7) 



(5.8) 



(5.9) 



(5.10) 



we can write 



S(s;tf) = .s* + -(2a s -s)* 2 + 0(* 3 ), 
S s (s;^) = -2a s *-fo s * 2 + 0(* 3 ). 



(5.11) 



Substituting the expansion for the metric, Eq. (3.2), and for the metric potentials, Eq. (4.1), we obtain to the 2.5PN 
order required for the equations of motion, 

1 + e (\v 2 G(l - QU a + GC(5 + 2s)U sa \ + e 2 H-v* + ^G 2 (l - () 2 U 2 
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+ \g{1 - Qv 2 U a - 4G(1 - (yv: + 1(5 + 2s)GCv 2 U sa - (5 + 2s)G 2 C(l - QU a U s 
+ |(15 + 18s + Aa s )G 2 t 2 U 2 sa + - ±JVi + i(5 + 2s)M>i 



+e 5/2 (W.5 + ^(5 + 2s)*i. 5 ) + 0(e 3 ) 



G = 

a 13 = 



p*v l 



l + e[^v 2 -G(l- QU a + GC(5 + 2s)U s 



Q(e 2 ) 



p*v l v j [1 + 0(e)] 
1 



* 2 



1 + e -w 2 - G(l - C)f/a + GC(5 + 2s)[/ SCT + 0(e 2 ) 



(5.12a) 

(5.12b) 
(5.12c) 
(5.12d) 



(1 - 2s) - e - 2s)v 2 + 3G(1 - 0(1 - 2s)U a - 3G( (l - 2s - ^a s J U, 
+ e 2 { - 1(1 - 2s)*; 4 + |G 2 (1 - C) 2 (l - 2s)U 2 ^G(l - 0(1 - 2s)v 2 U a 
+4G(1 - C)(l - 2syy; - |GC ( (l - 2s) - ^a s ) v 2 U sa - 9G 2 C(1 - ( 1 - 2s - ^a s ) [/„[/„ 



+ ^G 2 C 2 [l-2s 



+ £ 5/2 {^(l-2s-la s ) * 1 . 5 }+0(e 3 ) 



(5.12e) 



where U sa and are denned in terms of the cr-dcnsities. 

Substituting these formulas into the definitions of U a , U sa and the other potentials defined in terms of a, we can 
convert all potentials into new versions defined in terms of p*, plus PN corrections. For example, we find that the 
"Newtonian" potentials U a and U sa become 

U a = C/ + e|^$i-G(l-0$2 + 6GC$ 2s -GC$L 



4 ) + -G(l - C)S($i) + -G(l - Cmv 2 U) - 4G(1 - O^K^) - -G(l - £(*) 



-G 2 (l - C) 2 S($ 2 ) + ^G 2 (l - CfnU 2 ) + 9G&{v 2 U.) - ^G(X s (v 2 U s ) + l -GC, £.($?) 
-3GCS($^) + 3GCS(X S ) - ^GCS S (X S ) - G 2 ((l + 12Ai + 5C)£(*L) + G 2 C(2Ar + C)S a (*5.) 
+G 2 C(1 + 17C - 6A 1 )S([/ 2 ) - ^G 2 C(HC - 2\ 1 )Z S (U 2 ) - 6G 2 C(1 - QX(UU.) + G 2 ((l - QZ,(UU a ) 
+2G 2 ( 2 Z(a s U 2 ) - 6G 2 C(1 - C)S(*S) + G 2 C(1 - C)S 8 (*5) - 24G\ 2 £(£(a 8 Cg) + 4G 2 ( 2 £ s (£(a 5 C/ s )) 

f 4 (3) 1 ( 3 ) 1 ( 3 ) 1 

+£ 5/2 I _ j-kk {t) u __ X kk {t) {w -U,) + - 22 (t) (6^'C/ - ^C/ s - + ) - M s (i)(6£/ - U s )\ 
+Q(e 3 ), (5.13) 



- 3G(1 - C)$5 + 3GC$^ S - 4GCS(a s f/ s ) 



+e 2 --£ s (« 4 ) - -G(l - C)E fl (*i) - ^G(l - 0S s (« 2 C/) + 4G(1 - O^W - -G(l - C)£.(X) 



+G 2 (1 - 2 S,(* 2 ) + yG 2 (l - 2 S s (t/ 2 ) - \g& s (v 2 U s ) + 2G(Z(a s v 2 U s ) - ^GC^i) + 2G(Z(a s <£{) 
-2G(Z(a s X s ) + ^G(Z S (X S ) + 4G 2 C(2Ax + C)£(a.*2.) + G 2 C(1 - 4C - 6Ai)S a ($5.) - 4G 2 C(4C - A)£(a s [/ 2 ) 
+ ^G 2 C(2 + 7C - 6\i)Z s (U 2 ) - 4G 2 C 2 S(6 S (7 2 ) - 9G 2 ((1 - QV.{UU a ) + 12G 2 C(1 - CMa.UU,) 
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-12G 2 C 2 £ s (£(a s [/ s )) - 3G 2 C(1 - C)S 8 (*5) + 4 G 2 C(1 - C)S(a a *5) + ^G 2 ( 2 S(a s S(a s C/ s ))| 

f 1 (3) J (3) ^ 

+e 5 / 2 j -- If (t) (3(7 S - 4E(o 8 )) + - V a (t) (3x>U. - Ax>Y.{a a ) - 3X? + 4X(a 5 )' J ) - M s (t) (3U S - 4S(o 8 )) 

+0(e 3 ), (5.14) 
while the relevant PN potentials become 

<& lff = <!>! + e |^S(i; 4 ) - G(l - C)S(« 2 t/) + 6GOVig - G(Z s (v 2 U s )} + 0(e 2 ) , (5.15) 

<f 2 . = <f> 2 + e|^(« 2 C/) + ^S(<i> 1 )-G(l-OS([/ 2 )-G(l-OS(a> 2 ) 

+6GCS((7(7 S ) - G(Z S (UU S ) + 6GCS($ 2s ) - G<£($y j + 0(e 2 ) , (5.16) 

*2„ = ^+e{-is s (« 2 C/) + ^S s (cI. 1 )-3G(l-C)S s (C/ 2 )-G(l-0S s ($ 2 ) 

+3GCS S (C/C/ S ) - 4GC£(a s t/t/ s ) + 6GC£ s ($ 2s ) - GC£ a ($5,)} + 0(e 2 ) , (5.17) 

$ 2sff = 4> 2s + e|^E(« 2 t/ s ) - is($f) - G(l - C)S(t/C/ s ) - 3G(1 - C)E(*5) 

+6GCS((7 S 2 ) - GCS S ((7 S 2 ) + 3GCS(#5a) - 4GC£(£(a s t/ s ))j + 0(e 2 ) , (5.18) 

$1. = ^ s +e{-^(« 2 C/ s )-^ s (^)-3G(l-C)S s (?7?7 s )-3G(l-C)S 5 (^) 

+3GCS S (C/ 2 ) - 4G(Z(a s U 2 ) + 3GCE a (*L) - 4GC£ s (£(a s t/ s ))} + 0(e 2 ) , (5.19) 



X CT = X + e^-X(v 2 )-G{l-()X{U) + 6G(X{U s )-G(X s (U s )j+0{e 2 ), (5.20) 
^ Sff = X s + e^-±X s (v 2 ) - 3G(1 - ()X S (U) + 3GCX S (U S ) - 4G(X(a s U s )} + 0(e 2 ) , (5.21) 
V* = V l + e \\n^v 2 ) - G(l - ()Vi + GGCVl - G^siv'Us) J + 0(e 2 ) , (5.22) 



where all potentials are now defined in terms of the density p* , and including, where needed, the sensitivity factors 
s, a s and b s . In manipulating these expressions, we have made use of the identities, valid for any function /, 
E(s/) = [£(/) - S s (/)]/2 and E(a; 4 /) = x l Y>{f) - X'\f). The potentials U and C/ s will henceforth be given by 

= I (l-^(x-)),-( t ,x- )<iV 
/m |x-x'| 

In some cases we will use the same notation as before, to avoid a proliferation of hats, tildes or subscripts. We redefine 
the S, X and Y potentials by 

£(/) ee / P * (f ' X 2 /( *' x/) dV = P(4V/), (5.24a) 
Ja^ l x x I 

E'(/) ee / P * ft 'f K {^ X/) rfV = P(4 7 rp^7), (5.24b) 
Ja^ l x _ x I 

s ,. (/) ^ f ^( t »y p,J /( t ^) rf» a / = p ( 4 ff p' t ,V/) > (5.24c) 

./.M l x ~ x I 



(1 - 2s(x'))p*(t,x')/(i,x') , . , , 

i k JHy ' JJy ' J d 3 x' = P(4tt(1 - 2s)p*/) • 

X |x x I 



X(f) = f p*(t,x')/(i,x')|x-x'|dV, 

y(.f) = / P *(t,x')/(t,x')|x-x'| 3 dV 
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(5.24d) 
(5.24e) 
(5.24f) 



and their obvious counterparts X 4 , X s , F 4 , F 4 - 7 , F s , and so on. With this new convention, all the potentials 
defined in Eqs. (3.13) can be redefined appropriately. 



B. Equation of motion in terms of potentials 

Pulling together all the potentials expressed in terms of p*, inserting into the metric, Eq. (3.2), calculating the 
Christoffel symbols, we obtain from Eq. (2.5) the equation of motion 



where 



dv l /dt = a% + ea PN + e 3/2 a\ 5PN + e 2 a 2PJV + e 5/ V 2 5PN + 0(e 3 ) , 
a% = G(1-C)^+GC(1-2 S )C/' 4 , 



(5.25) 
(5.26) 



a PN = v 2 [G(l - Of/'* - G((l - 2s)U- 1 ] - 4G(1 - <>Vt/' j - v* 3G(1 - QU - G((l - 2s)U t 
-4G 2 (1 - Q 2 UU^ - 4G 2 C(1 - 0(1 - 2s)UU£ - 2G 2 C [Ai(l - 2s) + 2(s'} U S U? 
+8G(1 - Qv j V^ + 4G(1 - C)^ + ^G(l - ()X- 1 + ^G((l - 2s)X£ 

+\g{\ - o$i i - \cai - 2 S nv l - g 2 (i - o 2 ^ - g 2 c(i - 0(1 - ^)®r 

-G 2 ( [1 - C + (2Ai + 0(1 - 2s)] $5/ - 4G 2 C 2 (1 - 2s)^(a s Us) , 



1 ( 3 ) 

a i.5PAf = g(l _ ^S) II ■ 



(5.27) 
(5.28) 



a 2PN 



4G(1 - Qv i v i v k V i ' k + v 2 v l [G(l - QU - GQ{1 - 2s)U : 



+v l v 3 



+v j v k 



4G 2 (1 - *2 + 4G 2 C(1 - 0*L J - 2 G(1 - - 2G(1 - 0* 
2G(1 - C)$f' 1 - 4G(1 - ()$\ j ' k + 2G 2 (1 - CfPi^ - 4G 2 (1 - Q 2 P l 2 j ' k 



+2G 2 C(1 - 0^'f - 4G 2 C(1 - 0^ 



2 s 



+« 2 



--G(l - 0*i' + -GC(1 - 2s)$f 4 - G 2 (l - C) 2 * 2 l + G 2 C(1 - 0(1 - 2s)$^ 
-G 2 C [1 - C - (2Ai + 0(1 - 2s)] $£/ + 2G 2 C [Ai(l - 2s) + 2(s'} U.U? 



1 



1 



+4G 2 C 2 (1 - 2s)£- 4 (a s [/ s ) + -G(l - ()X' 1 - -G((l - 2s)X_ 
3G 2 (1 - C) 2 * 2 - G 2 C(1 - 0(1 - 2s)^ + G 2 C [3(1 - - (2Ai + 0(1 - 2s)] <&L 



1 



1 



-4G 2 C 2 (1 - 2s)£(a s £/ s ) - 2G 2 C [Ai(l - 2s) + 2C«'] U S U S - -G(l - 0*1 - ^G((l - 2s)$? 
~G{1 -C) X +2^C(l - 2s) X s +4G 2 (1 - C) 2 V k U< k + 4G 2 C(1 - 0(1 - ^)V k V k 
+v j [8G 2 (1 - () 2 v¥' j] + 8G 2 C(1 - O^'V 1 ^) - !6G 2 (1 - Q 2 <S>f i] + AG ( l ~ C)* M 
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+32G 2 (1 - () 2 G [ j' j] - 8G 2 ((1 - ()P(U S U^ ] - 16G 2 (1 - Q 2 UV^ - 4G(1 - 0£ ,[i (vV) 
+8G 2 (1 - () 2 V l U' J + 8G 2 C(1 - C)(l - 2s)VW> 1 - 4G(1 - 0*? - 4G 2 (1 - () 2 P? - 4G 2 C(1 - QP^ 

1 (4) i (4) (3) o -, 

+ -G(1 - C) +^GC(1 - 2s) Yf +2G(1 - C) X* +-G(1 - 0*i* - -GC(1 - 2s)l<> 2 ) 
+2G(1 - ()±(v*v 2 ) + 7 -G(l C)S^(« 4 ) - ^GC(1 - 2s)E>> 4 ) + ^G 2 (l - () 2 ^(v 2 U) 
-1g 2 C [3(1 - C) - (2Ai + C)(l - 2a)] ^(v 2 U s ) ^G 2 ((l 0(1 2s)£>> 2 [/) 

+2G 2 C 2 (1 - 2s)^\v 2 a s U s ) - 4G 2 (1 - C) 2 S ,i (w 3 V J ') + 4G 2 ((1 - 0(1 - 2«)Ej,VV J ') - ^G 2 (l - O^'^i) 

-^G 2 C(1 - 0(1 - 2 S )S^$i) + 2G 2 C 2 (1 - 2s)£> s $ s 1 ) + 1g 2 C [1 - C + (2Ax + 0(1 - 2s)] £;*(*;) 

-6G 2 (1 - 2 ^*i + 2G 2 C(1 - 0(1 - 2«)I7$; ,i + G 2 C [Ai(l - 2s) + 2(s'] ?7 s $f l 
-2G 2 (1 - 2 $i^' 4 - 2G 2 C(1 - 0(1 ~ 2s)$it/' 1 + G 2 C [Ai(l - 2s) + 2(s'] &\U£ 
-4G 2 (1 - () 2 $\ j U' j - 4G 2 C(1 - 0(1 - 2s)$\ j U s j + 8G 2 (1 - QfVW^ 
+4G 2 (1 - () 2 V l U - 4G 2 C(1 - 0(1 - 2s)V l U s 

-2G 2 (1 - Q 2 UX^ - 2G 2 C(1 - 0(1 - 2s)UX? - G 2 ( [Ai(l - 2s) + 2(s'] U 8 X£ 
-2G 2 (1 - CfXU' 1 - 2G 2 C(1 - 0(1 - l8)XU'j - G 2 C [Ai(l - 2s) + 2(s'] A J ;' 

-8G 2 (1 - QfUV' - ^G 2 (l - 2 ^\X) - iG 2 C(l - 0(1 - 2s)^:(X) 
-\g\ [1 - C + (2Ai + 0(1 - 2s)] £?(*.) - 2G 2 ( 2 (1 - 2s)V i {a s X s ) 
- l -G 2 {l Q 2 Xi 2G 2 C 2 (1 - 2s)X> i (a s Us) ^((l - 0(1 - 2s)X-'(C/) 
-1g 2 C [1 - C + (2Ax + 0(1 - 2s)] X;;(U S ) + 4G 2 (1 - Q 2 \q - 4G 2 ((l - QVtf 

-8G 2 (1 - 2 $ 2 - 6G 2 (1 - () 2 G'I + 2G 2 C(1 - OG'L - 4G 2 (1 - C) 2 G^ - 4G 2 ((1 - 0(1 - 2s)G 2 4 s 
+8G 2 (1 - Q 2 G'Z + 8G 2 C(1 - 0(1 - 2s)G^ s + 8G 2 (1 - Q 2 G'l 

-4G 2 (1 - () 2 Gg - 4G 2 C(1 - 0(1 - 2s)G^ s + 16G 2 (1 - 2 G^ - 4G 2 ((1 - QP(U a U?) 
+4G 3 (1 - Q 3 U*i + 4G 3 C(1 - [1 - C + (2Ai + 0(1 - 2s)] U$ s 2 / + 4G 3 ((1 - 2 (1 - 2s)U^ 1 
+2G 3 C(1 - [Ai(l - 2s) + 2(s'] U&tf + 2G 3 ((2\ 1 + () [Ai(l - 2s) + 2(s'] [/,$»/ 
+16G 3 C 2 (1 - 0(1 - 2s)UY,' l {a s U s ) + 8G 3 ( 2 [Ai(l - 2s) + 2(s'] £/ a E> i (a a I7 a ) 
+4G 3 (1 - Cf^U' 1 + 4G 3 C(1 - Cf^LU' 1 + 2G 3 C(1 - [Ai(l - 2s) + 2(s'] &' 2 U$ 
+2G 3 C{2C(1 - 0(1 - 2s) + (2Ai + 0[Ai(l - 2s) + 2Cs']} + 4G 3 C(1 - C) 2 (l - 2*)$ 2 I/- i 

+8G 3 C 2 [Ai(l - 2s) + 2(s'] Y.(a a U a )U<j 

+8G 3 (1 - () 3 U 2 U> 1 + 8G 3 C(1 - 2 (1 - 2s)U 2 U£ + 8G 3 ((1 - [Ai(l - 2s) + 2(s'] UUJJ'j 

+G 3 C [(8A 2 - 2(Ai - 2A 2 )(1 - 2s) + 12AiC«' - 4C 2 s"] U 2 U? 

-G 3 (l - 3 £"(* 2 ) - G 3 C(1 - O^'^L) - G 3 C(1 - 2 (1 - 2s)Ei l ($ 2 ) 

+G 3 C {(2Ai + [1 - C + (2A a + 0(1 - 2s)] - C(l - 0(1 - 2s)} E^L) 

+G 3 C(1 - [1 - C + (2Ai + 0(1 - 2s)] E;($5) + 4G 3 C 2 (1 - 0(1 - 2s)E- i (a s ^) 

+4G 3 C 2 (2A! + 0(1 - 2s)E' i (a a $5 a ) 

+16G 3 C 3 (1 - 2s)E> i (a s E(a s t/ s )) + 4G 3 ( 2 [1 - C + (2Ai + 0(1 - 2s)] E' 4 (E(a s C/ s )) 
+ ^G 3 (1 - () 3 ^{U 2 ) + ^G 3 C(1 - 2 (1 - 2s)^ s l (U 2 ) + G 3 C(1 - 2 ^(^ 2 ) 

+ ^G 3 C {(2Ai + 0(1 - + (1 - 2s) [C(l - + C(2Ai + 1) + 16A 2 - 4A 2 ] } £^(t/ 2 ) 

+G 3 C(1 - [1 - C + (2Ai + 0(1 - 2s)] ^(U S U) + 2G 3 C 2 [1 - C + 6Ai(l - 2s)] £'> S C/ S 2 ) 
-4G 3 C 3 (1 - 2s)E- i (6 s f/ 2 ) + 4G 3 C 2 (1 - 0(1 - 2s)E' i (a s U s U) 
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l 2.5 PN 



-4G 3 (1 - () 3 P?U> j - 4G 3 C(1 - CfPllU' 3 - 4G 3 C(1 - C) 2 (l - tyPgU? - 4G 3 C 2 (1 - 0(1 - tyPgU? 
-4G 3 (1 - CfH' 1 - 4G 3 C(1 - CfH^ - 4G 3 C(1 - C) 2 (l - 2s)ff s ■* - 4G 3 C 2 (1 - 0(1 - 2s)i? s s,< , (5.29) 

o /(5) 1 (5) \ (4) (3) 

- x i to gij x kk + 2v j T* +2 [G(l - t)U- j + G((l - 2s)U'J ] V 3 

4 (3) (3) 

+ - [G(l - ()U' 1 + GC(1 - 2s)U£] I kk - [G(l - ()X^ k + GC(1 - 2s)X^ k ] V k 

2 (5) 2 (4) 1 / (5) 1 (5) \ 1 / 1 \ (5) 

-- T» +-e™ J« --(1 - 2s)^' f +-«*« Xf J + -(1 - 2s) (sV + -r 2 ^j Xj 

1 (5) i (3) 4 (3) 

+-(1 -2s) T» -_ w 2(i_2s)Xi--G(l-C)(l-2s)C/li 

1 . / . W W .. \ 1 . (3) 

+ - 2s) 2x ] If - I kk -6M S - -(1 - 2s)x l M s 



i r i / (3) (3) 

--G|[l - C + (4Ai + 0(1 - 2s) + 4Cs'] U? + 4C(1 - 2s)X>\a s )\ f 2x j V s - if -6M S 
~ 3 G { [1 ~ C + (4Al + C)(1 ~ 2S) + 4C ' S ' ] Us + 4C(1 ~ 2s ^ a ^ } V * 



1 



(3) 



+-G|[l-C+(2A 1+ 0(l-2s)]X^+4C(l-2s)X^(a s )|Xf . (5.30) 

We next turn to the problem of expressing these equations explicitly in terms of positions and velocities of each body 
in a two-body system. 



VI. EQUATIONS OF MOTION FOR TWO COMPACT BODIES 

We now wish to calculate the equation of motion for a member of a compact binary system. To do this, we integrate 
p*dv l /dt over body 1, and substitute Eq. (5.25) and then Eqs. (5.26) - (5.30). We follow closely the methods already 
detailed in [26] (hereafter referred to as PWII) for evaluating the integrals of the various potentials, and so we will 
not repeat those details here. Readers should consult Sec. Ill and Appendices B, C, and D of PWII for details. In 
structural terms almost all of the potentials that appear in the 2PN terms in scalar-tensor theory also appear in 
general relativity, apart from the differences in the types of densities that generate the potentials, for example U s vs. 
U, X s vs. X, vs. $ 2 , and so on. The only 2PN term that does not appear in GR involves the potential P(U S U£), 
but this can be evaluated using the methods described in PWII. 

Similarly, at 2.5PN order most of the moments that appear here also appear in GR, only a few, notably the 
scalar monopole and dipole moments M s and I\ are new. Particularly new is the appearance of a 1.5PN order term 
generated by the scalar dipole moment; this, of course, is the radiation-reaction counterpart of the well-known dipole 
gravitational radiation prediction of scalar-tensor theories. 

A. Conservative 1PN and 2PN terms 

We begin with the conservative Newtonian, 1PN and 2PN terms. The results are, at Newtonian and 1PN orders. 

<4 ( pjv) ~ n +— ^n l |-(l + 7K-(2 + 7)(« 2 -2 Vl .v 2 ) + -(v 2 .n) 2 

+ [4 + 27 + 2/3i] + [5 + 2 7 + 2/3 2 ] — — - \ 

+ Gam 1{vi _ [(4 + 2 _ )vi . n _ (3 + 2 ^ )va . n] > 
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4(PiV) = (1-2), (6.1) 

where r = |xi — X2I, n = (xi — X2)/r, and where the parameters a, 7, and /3a are defined in Table I. Note that under 
the interchange (1^2), n4 n. At 2PN order, we find 

°W) = ^ 2 ^n 4 {-(2 + 7)h 4 -2 U2 2 (v 1 .v 2 ) + (v 1 .v 2 ) 2 +3(v 2 .n) 2 (v 1 .v 2 )] 
+ + 7> 2 (v 2 • n) 2 + I (3 + 7) «2 2 (v 2 • n) 2 - y (v 2 • n) 4 

+ G^m2 / (2 + _ } ^ _ _ ^ _ 2 + I (( 2 + 7) 2 + 4( f 2 ) [(vi • n) 2 - 2(vi • n)(v 2 • n)] 
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,2 

^ (57 + 447 + 97 2 + 16(3 + 7)^2 + 4Ji - 8x2) 
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1 G 2 a 2 m 2 



"2 
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+ ^K-^2(2 + 7) 
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+ (1 + 7)« 2 (v 2 • n) - (5 + 37)w 2 (v 2 • n) + |(3 + 27)(v 2 • n) 3 

+ ( (55 + 407 + 27 2 - 16^2 + 8J1) v 2 • n - (63 + 407 + 2^ 2 - 8/3 2 + 85i) vi • n 

( ((2 + 7) 2 + 45 2 ) Vl • n + (4 - 7 2 + 4/3! - 4S 2 ) v 2 • J , 



1 Gam 2 



2 r 

4(2PiV) = (1-2), (6.2) 

where 5a and xa are defined in Table I. 

It is straightforward to show that these equations of motion can be derived from a two-body Lagrangian, given by 

1 2 1 4 1 6 \ i 1 Gamim-2, 
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+ | + 2 ^ _ I (? + 4 ^ )V1 ' V2 _ I (V1 ' n)(v2 ' n) _ 2 (1 + 2/?2) 

+ (I(7 + 47) K - , 2 (v 2 • n) 2 ] - (2 + 7K( Vl • v 2 ) + i( Vl • v 2 ) 2 
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2 1 , -2-4 
1+ g7+g7 +2/3 2 + -5i--X2 



+ - 



1 G 2 a 2 miiri2 



[19 + 87 + 8^1 + 8^2 - 327-^1^2] 
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— -Gam 1 ra 2 
8 

+ (1-2). 



2(7 + 47)ai • v 2 (v 2 ■ n) + n ■ ai(v 2 • n) 2 - (7 + 47)11 • ait> 2 



(6.3) 



As in general relativity, the Lagrangian contains acceleration-dependent terms at 2PN order, and thus the Eulcr- 
Lagrange equations are (d 2 /dt 2 )(5L/5a t ) — (d/dt)(8L/5v l ) + SL/5x l = 0. The equations of motion (absent radiation- 
reaction terms) admit the usual conserved quantities. The energy is given to 2PN order by 
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while the total momentum is given by 
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B. Radiation- reaction terms 

At 1.5PN order, the leading dipole radiation reaction term is given by 

1 W 
a l(1.5PJV) = g(l~ 2s l)^s; 

i (3) 

4(i.5Piv) = 3(1-2*2)1; . (6.6) 

Because we will be working to 2.5PN order, the scalar dipole moment I* must be evaluated to post-Newtonian order, 
and when time derivatives of that moment generate an acceleration, the post-Newtonian equations of motion must 
be inserted. Explicit two-body expressions for 2] and the other moments needed for the radiation-reaction terms are 
provided in an Appendix. In addition to evaluating the direct 2.5PN terms from Eq. (5.30) for two bodies, we must 
include the 1.5PN contributions to the accelerations that occur in the 1PN terms V 1 , X' 1 and Xg that appear in Eq. 
(5.27). 

At 2.5PN order, the final two-body expressions take the form 

o /(5) I (5) \ (4) 1 r (3) r (3) „ (5) 2 (4) 
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+-v\(l - 2*0 2zj If - If -6M S - 3(1 - 2*1)4 M fl - g«i(l - 2*i) T B 



1 Gumi 
6 r 2 



f 1 / • (3) (3) \ 

'1-2,- 47- 1 [(1 - 2si)& + (1 - 2« 2 ) j 8 2 ] I I 2^ If - If -6M S J 

-^nV(l - 2 S2 )(1 - 8^/7) S -^(1 - 2 fll )(l - 8/V7) S§ 

+|^(«i- a )(7 + 47)S, 
4(2. 5 Piv) = (1-2). (6.7) 

We shall defer calculating the moments and their time derivatives explicitly until the next subsection, where we obtain 
the relative equation of motion. 



C. Relative equation of motion 



We now wish to find the equation of motion for the relative separation x = xi — x 2 , through 2.5PN order. We take 
the PN contributions to the equation of motion for body 1 and body 2 and calculate cPn/dt 2 = &i — a 2 . We must 
then express the individual velocities vi and v 2 that appear in post-Newtonian terms in terms of v = vi — v 2 . Since 
velocity-dependent terms show up at 1PN order, we need to find the transformation from vi and v 2 to v to 1.5PN 
order so as to keep all corrections through 2.5PN order. To do this we make use of the momentum conservation law 
(6.5). But because of the contributions of dipole radiation reaction at 1.5PN order, the momentum is not strictly 
conserved because of the recoil of the system in response to the radiation of linear momentum at dipole order. By 
combining Eqs. (6.5) and (6.6), it is straightforward to show that the following quantity is constant through 1.5PN 
order: 

miv\ (l + \v 2 ^j - \ Gam ^ [„j + n i( Vl . n) ] + I mi (i _ 2si )±l + (1 ^ 2) = & . (6.8) 
Setting C % = and combining this with the definition of v, we find that 



25 



v\ = -—!;* + **, (6.9) 
m 



where 



1 , X ( o Gam\ . Gam . -1 2 > _ , _ , _ . / Gam \ 2 ,• _ . 9 . ,„ „ „. 

-r?V U 2 -— — J «' - -j^-fa;' - - C ^ + +^_) f — 1 »* + 0(e 2 ) , (6.10) 



where m and are the total mass and reduced mass ratio, -0 = Sm/m = (mi — m 2 )/m, and 

S_ = -a- 1/2 ( Sl - S2 ), 

S+ = a- 1/2 (l- Sl -s 2 ). (6.11) 

We also need to evaluate the multipole moments that appear in the radiation-reaction terms to the appropriate 
order, and then calculate their time derivatives, inserting the equations of motion to the appropriate order as re- 
quired. Explicit formulae for the moments are displayed in Appendix A. Combining all the various PN contributions 
consistently, we arrive finally at the relative equation of motion through 2.5PN order, as given in Eqs. (1.4) and (1.5). 
Here we display the 2.5PN coefficients: 

2 Gam 2 

A-2.5PN = a\v + a 2 ha 3 f , 

r 

„ 9 Gam , .2 .„ „ „. 

B 2 .5PN = hv 2 + b 2 + b 3 r 2 , (6.12) 



where 



5 15-5 15 

ai = 3 - -7 + y/3+ + gC5 2 (9 + 4 7 - 2rf) + , (6.13a) 

« 2 = y + f 7 - f ^ + - t 135 + ^ + 8 ^ + 32 ^ + ] + ( S - h+ - S+ ^ 

-§C*S- (* - f ^- + 16 5 ^;^- ) - 40 C (&^) 2 , (6,3b) 
25 

« 3 = y & - ( X - 2r ?) - - C^-5+] , (6.13c) 

&i = l-^7 + ^/3+-^C5 2 (7 + 47-2r / ) + ^C^_5 + , (6.13d) 

fe 2 = 3 + ^7 - §0+ - ^CSl [23 + 87 - 8, + 8/3+] + y tf- (^y^ 

I + ~ 3 ^" Y 7 J ' (6.13c) 

b 3 = § [67 + (S 2 _ (13 + 87 + 2rj) - 12P+ - 3^SS + ] . (6.13f) 



D. Energy loss rate 

We now wish to evaluate the rate of energy loss that is induced by the radiation-reaction terms in the equations 
of motion. Because those equations of motion contain both 1.5PN as well as 2.5PN contributions, we will have not 
only the normal "quadrupole" order contributions to the energy loss rate analogous to those that appear in general 
relativity, but also dipole contributions that are in principle larger by a factor of 1/v 2 . Since the conventional "counter" 
for keeping track of contributions to the waveform and energy flux in the wave- zone denotes the GR quadrupole terms 
as "Newtonian" or OPN order, the dipole terms will, by this reckoning, be of — 1PN order. 

To evaluate the energy loss correctly through "Newtonian" order, we first express the conserved energy in relative 
coordinates to 1PN order. Using the transformations (6.9) and (6.10) to 1PN order, we obtain 

„ 1 9 Gam 3 „ . . 

E = -iiv 2 - ix—- + - 3t7> 4 
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1 Gam r .„ , 9 . 9n 1 / Gam\ 2 , „ - „ , - . . _ 

+ — [(3 + 27 + r,) V 2 + r;r 2 ] +-^f — — J (1 + 2/3+ - 2^/3_) . (6.14) 

We then calculate dE/dt, inserting the 1.5PN and 2.5PN acceleration terms into the leading term v • a, and inserting 
only the 1.5PN terms wherever accelerations occur in the time derivative of the 1PN terms. 

Beginning with the leading term, and expressing the 1.5PN acceleration in the form &1.5PN = (D/r 3 )(3rn — v), 
where D = ArjC{Gam) 2 S 2 _/3, we find for the -1PN term (dE/dt)- 1PN = n(D/r 3 )(3r 2 - v 2 ). This can be simplified 
by exploiting the identity 

d ( f \ v 2 - 3f 2 + x • a 

-?. • ( 6 - 15 ) 



dt \ r 2 

Thus (v 2 — 3r 2 )/r 3 can be written as the total time derivative of a quantity that can be absorbed as a 1.5PN correction 
to the definition of E, leaving (dE/dt)-\PN = fi{D/r 3 )(x ■ a). Inserting the Newtonian acceleration for a, we obtain 

{dE,dt UpN = -\^(^S 2 _. (6.16) 

This is in agreement with earlier calculations of the energy flux due to dipole gravitational radiation [36, 37]. 

However, since we are working to Newtonian order in the energy loss, we also need to include the 1PN contributions 
to the acceleration that appears in Eq. (6.15), yielding a contribution given by /iD(Gam/r 4 )(ApN + r 2 Bpjy), where 
Apn and Bpn are given by Eqs. (1.5a). We then combine this with the other Newtonian order terms generated from 
dE/dt, leading to an expression of the general form 



dE 8 jiiq ( Gam\' 
~db ~ 157" V r ) 



Gam 2 Gam , 2 2 . 2 f Gam\ 2 4 . 4 
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Pi v +P2 r +P3V r +p 4 +P5V +p&r 

r r 

We now use an identity derived from the Newtonian equations of motion, 

pv 4 - pv 2 ^^ -(p+ q)v 2 r 2 - 2s^^r 2 ) . (6.18) 



d f v 2s r p \ _ v 2s 2 r p 1 / 4 2 Gam 2 ., Gam , 2 



dt\ri J ri +1 

This is applicable at this PN order provided that the integers s and p are non-negative, q > 2 and 2s + p + 2q = 7. 
Using the three possible cases (s, p, q) = (1, 1, 2), (0, 3, 2), (0, 1, 3), we can freely manipulate the values of three of 
the six coefficients pi in Eq. (6.17). The idea is to combine terms on the right-hand-side of Eq. (6.17) into a total 
time derivative, to move that to the left-hand-side and then to absorb it into a meaningless redefinition of E (see for 
example, [50, 51] for discussion). Thus one can easily arrange for p^, p§ and pe to vanish. It then turns out that the 
coefficient p^ of the term proportional to v 2 r 2 is proportional to the combination of the 2.5PN equation-of-motion 
coefficients 5a\ + 3a$ — 156i — 563. An inspection of Eqs. (6.13) reveals that this combination miraculously vanishes. 
Pulling everything together, we obtain the final expression for the energy loss rate, 

(dE/dt) 0PN = -A^ fGamV _ ^ . 2) 

15 r \ r / 
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+40 C ^ - 30C*5- ( S+h + S -' P -) - 120C ( S+h+ - S -' P -) 2 ■ 
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1 

(6.20) 



These results are in complete agreement with the total energy flux to —1PN and 0PN orders, as calculated by 
Damour and Esposito-Farese [35]. (We are grateful to Michael Horbatsch for his invaluable help in verifying this 
agreement.) 
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TABLE II: Dictionary of parameters used in the equations of motion. DEF refers to Ref. [20, 35]; TEGP refers to Sec. 11.3 of 
Ref. [37] ; PPN refers to the parametrized post-Newtonian limit of weakly gravitating bodies 



This paper 


DEF 


TEGP 


PPN limit 


Ga 


G12 


Q12 


1 


7 


712 


|(Bl2/ffl2 - 1) 


7-1 


fa 


022 


!(2Wei 2 2-i) 


/3-1 


fa 


#1 


i(2?21l/ai 2 2-l) 




Si 








8 2 








Xi 


1,1 

~ 4^222 






X2 


~4 e lll 






T^ifa 


— |Cl212 







VII. DISCUSSION 

We have used the DIRE approach based on post-Minkowskian theory to derive the explicit equations of motion in 
a general class of massless scalar-tensor theories of gravity for compact binary systems through 2.5PN order. Here we 
discuss the results, and compare our work with related work on scalar-tensor gravity and equations of motion. 

A. General remarks and comparison with other results 

We begin by noting that, not surprisingly, the expressions are considerably more complicated than the corresponding 
general relativistic expressions. Given that the results depend on the function u(<j>) and its first and second derivatives, 
on the masses of each body, and on the sensitivities of each body and their derivatives, it is somewhat remarkable 
that the final equations of motion depend on a rather small number of parameters, as shown in the right-hand column 
of Table I. The parameter a combines with G to yield an effective two-body Newtonian coupling constant. It is 
not a universal constant, as it depends symmetrically on the sensitivities of each body. The parameter 7 and the 
body-dependent parameter 0a govern the post-Newtonian corrections, while the body-dependent parameters 5a and 
Xa govern the 2PN corrections. In the radiation-reaction terms, the sensitivities sa occur explicitly along with 7 and 
Pa- 

The relative simplicity of the parameters at 1PN and 2PN orders has been noted before. Damour and Esposito- 
Farese [20, 35] (DEF hereafter) studied a class of multi-scalar-tensor theories, but worked in the Einstein representa- 
tion, where the gravitational action was pure general relativity, augmented by a free action for the scalar fields. This 
is a non-metric representation of the theory, since the scalar field(s) couple to normal matter via a function A(<p) 
(here we will focus on a single scalar field). For a compact body with mass m(<p) (using the Eardley ansatz), the 
effective matter action depends on the product A(ip)fh(ip). The scalar field 4> of our Jordan representation is given 
by <fi = A((p)~ 2 , and 3 + 2uj(<j>) = (d In A/ dip) ~ 2 . Using a diagrammatic approach, DEF showed that the important 
quantities involved derivatives of A(ip)m(ip) with respect to <p, and consequently (in our language) uj and sa and their 
derivatives always combined in specific ways, leading to relatively few parameters. Table II gives a dictionary that 
translates from our parameters to those of DEF for the case of two bodies. Interestingly, our parameters 5a do not 
appear in DEF's list, so far as we could tell. 

In the 1PN limit, Will [37] wrote down a general A^-body Lagrangian for compact self-gravitating bodies that could 
span a wide class of metric theories of gravity that embody post-Galilean invariance (so-called "semi-conservative" 
theories of gravity), and that have no "Whitehead" potential in the post- Newtonian limit. Comparing our Lagrangian 
of scalar-tensor theory with the 2-body limit of Eq. (11.62) of [37], we can translate between our parameters and the 
coefficients Gab, Bab, and T> a b c of [37], as shown in Table II. 

The factor 1 — 2s a appears throughout these equations. This quantity is often called the "scalar charge" of the 
object. From the point of view of the Einstein representation of scalar-tensor theory, it is easy to see how this factor 
arises. The scalar field appears in the gravitational part of the action only in a kinetic term g^ip^ip^ (we assume 
that there is no potential V(<p)). It docs not couple to gravity other than via the metric in the kinetic term. The 
effective matter action for a compact body depends on the product A(ip)M(ip). Varying this product with respect to 
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ip yields the quantity 

a / \ , r , \ ( d\wA d\wM d\n<h\ din A .„ „ . „ . . 

^)M(„) ^_ + — ^ j ^ = - 2*)^ , (7.1) 

where we used the fact that ln</> = — 2 In A(<p). Thus the factor 1 — 2s and its derivatives naturally control the source 
of the scalar field, as can be seen clearly in Eq. (5.12e). Defining a scalar charge for body A in a two-body system by 

q A = a- 1 ' 2 {l-2s A ), (7.2) 

we see that the quantities S± are given by 

S+ = ^(qi + 92) , 

S- = -92). (7.3) 

The scalar charge, or sensitivity of a given body depends on its internal structure. For weakly gravitating bodies, 
s rts —Sl/M <C 1, where f2 = — (1/2)G J p*p'*\x — x'\~ 1 d 3 xd 3 x' is the Newtonian self-gravitational binding energy . 
For neutron stars, values of the sensitivities range from 0.1 to 0.3, depending on the mass and equation of state of 
the body [52, 53] and can vary dramatically, depending on the specific form of uj{4>) [35]. 



B. Weakly self-gravitating systems 

In the post-Newtonian limit with weakly self-gravitating systems, the sensitivities are themselves of order e. If 
one is working purely at 1PN order, then the effects of sensitivities in the 1PN terms of Eq. (6.1) will be of 2PN 
order. So the only effect of the bodies' sensitivities in this case will come from the coefficient a in the Newtonian 
term. Consider a specific example: body 1 with sensitivity s\ resides in the field of body 2, with sensitivity zero. The 
acceleration of body 1 is then given by 

ai = -^V(l-2C*i), (7.4) 

and thus the body's Newtonian acceleration will depend on its internal structure, a violation of the Strong Equivalence 
Principle, commonly known as the Nordtvedt effect. In the PPN framework [37], the Nordtvedt effect is normally 
expressed in terms of O. Alternatively, since M w mo + 0, we have that O/M = dlnM/dlnG. Taking into account 
Eq. (4.6), we can connect the sensitivity s to by 

where 

~ (4 + 2 Wo ) 2 (3 + 2 Wo ) 1 ■ ' 

is the parameter defined in TEGP (see Eqs. (5.36) and (5.38)) such that the PPN parameter (3 = 1 + A in scalar-tensor 
theory (note the relationship between (j) and G, which is set equal to unity in TEGP). We also have that 7 = 1 — 2£. 
We can then express the acceleration of body 1 as 



Gm 2 
»i = 5-n 



l + (-^ + 4A^ 



(jJq J mi 



(7.7) 



The coefficient in front of J7i/mi is precisely 4/3 — 7 — 3, as in the standard PPN framework. 

In the 1PN terms in Eq. (6.1), for weakly self-gravitating systems, it is easy to see from Table I that in the limit 
Si — > 0, a — > 1, the parameters 7 and /3j tend to the PPN parameters 7 — 1 and /3 — 1, respectively, as shown in Table 
II, and thus our equations of motion at 1PN order agree with the standard ones for "point" masses in scalar-tensor 
theory. 

The radiation-reaction results can also be compared with existing work. The —1PN energy loss due to dipolc 
gravitational radiation reaction, Eq. (6.16) is in complete agreement with calculations of the dipole energy flux [36, 
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37, 54]. In comparing Eq. (6.16) with Eqs. (10.84) and (10.136) of [37], the additional factor of [1 +4A(2 + cj )] 2 arises 
from the relation (7.5) between s and Q./M. 

For weakly self-gravitating bodies, the Newtonian-order energy loss simplifies by virtue of setting all sensitivities 
equal to zero. In this case, with a = 1, 7 = — 2(, f3 + — (3 — 1 = A, /3_ = 0, S- = 0, and S + = 1, wc obtain 



m = 12 

K2 = 11- 

= 11- 



2 + cjo 
45 



C - 40A - 30A7C 



45 



8 + 4uj 

These agree completely with Eq. (10.136) of [37]. 



1 8 /2A 

1+ Ht 



2A 

T 



(7.8) 



C. Binary black holes 



Roger Penrose was probably the first to conjecture, in a talk at the 1970 Fifth Texas Symposium, that black holes in 
Brans-Dicke theory are identical to their GR counterparts [55] . Motivated by this remark, Thorne and Dykla showed 
that during gravitational collapse to form a black hole, the Brans-Dicke scalar field is radiated away, in accord with 
Price's theorem, leaving only its constant asymptotic value, and a GR black hole [55]. Hawking [16] proved on general 
grounds that stationary, asymptotically flat black holes in vacuum in BD are the black holes of GR. The basic idea 
is that black holes in vacuum with non-singular event horizons cannot support scalar "hair" . Hawking's theorem was 
extended to the class of f(R) theories that can be transformed into generalized scalar-tensor theories by Sotiriou and 
Faraoni [56]. 

For a stationary single body, it is clear from Eq. (5.12e) that, if s = 1/2 and all its derivatives vanish, the only 
solution for the scalar field is <fi = <po , and hence the equations reduce to those of general relativity. In the Einstein 
representation, this corresponds to A(ip)M((p) — constant, so that the scalar field decouples from any source, and 
thus must be either constant or singular. Consequently, stationary black holes are characterized by s — 1/2. 

Another way to see this is to note that, because all information about the matter that formed the black hole has 
vanished behind the event horizon, the only scale on which the mass of the hole can depend is the Planck scale, and 
thus M cx Mpi anck cx G~ 1/2 cx <t) 1/2 . Hence s = 1/2. 

If Sa = 1/2 for each black hole in a binary system, then, as we discussed in the introduction, all the parameters 
7, Pa, 5a, Xa, and S± vanish identically, and a = 1 — (. But since a appears only in the combination with GarriA, 
a simple rescaling of each mass puts all equations into complete agreement with those of general relativity, through 
2.5PN order. 

But is sa = 1/2 really true for binary black holes? If the orbital timescale is long compared to the dynamical 
(quasinormal mode) timescale of each black hole, then it is plausible to assume that Hawking's theorem holds for 
each black hole, at least up to some PN order. On the other hand, one could imagine a situation where each hole 
is distorted by the tidal forces from the companion hole, or where gravitational radiation flowing across the event 
horizons disrupts the stationarity needed for Hawking's theorem. In PN language, these kinds of effects are known 
to be of an order higher than the 2.5PN order achieved in this paper, so perhaps some non-GR effects might emerge 
at sufficiently high PN order. Can a perturbation of the scalar field be supported sufficiently by strong gravity or 
by time varying fields to make any difference? Or, without matter to support it, does any scalar perturbation get 
radiated away on a quasinormal-mode timescale, which is short compared to the orbital timescale, except during the 
merger of the two black holes? Preliminary evidence from numerical relativity supports the latter scenario: Healy et 
al. [57] introduced a very large Brans-Dicke type scalar field into the initial data of a binary black hole merger and 
found that, while the field affected the inspiral while it lasted, it was radiated away rather quickly, although it was 
not possible from the numerical data to fully quantify this. 

It should be pointed out that there are ways to induce scalar hair on a black hole. One is to introduce a potential 
V{4>), which, depending on its form, can help to support a non-trivial scalar field outside a black hole. Another 
is to introduce matter. A companion neutron star is an obvious choice, and such a binary system in scalar-tensor 
theory is clearly different from its general relativistic counterpart (see the next subsection). Another possibility is a 
distribution of cosmological matter that can support a time- varying scalar field at infinity. This possibility has been 
called "Jacobson's miracle hair-growth formula" for black holes, based on work by Jacobson [58, 59]. Whether it is 
possible to incorporate such ideas into our approach is a subject for future work. 

These considerations motivate us to formulate a conjecture along the following lines: Consider a scalar-tensor theory 
of gravity with no potential for the scalar field, and consider two black holes with non-singular event horizons in a 
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vacuum (no normal matter), asymptotically flat spacetime with <f> at spatial infinity constant in time. Following an 
initial transient period short compared to the orbital period, the orbital evolution and gravitational radiation from 
the binary system are identical to those predicted by GR, after a mass rescaling, independent of the initial scalar field 
configuration. Aspects of this conjecture could be addressed by numerical simulations that extend the work of [57]. 
It may also be possible to address it partially by generalizing Hawking's theorem to a situation that is not strictly 
stationary, but yet still retains some symmetry, such as a helical Killing vector. This will be the subject of future 
work. 



D. Black- hole neutron-star systems 



Finally, we note the unusual circumstance that, if only one of the members of the binary system, say body 2, is a 
black hole, with s 2 = 1/2, then a = 1 — £, 7 = 0a = 0, and hence, through 1PN order, the motion is again identical 
to that in general relativity. This result is actually implicit in the post-Newtonian equations of motion for compact 
binaries in Brans-Dicke theory displayed in Eq. (11.91) of [37], but was never stated explicitly there. 

At 1.5PN order, dipole radiation reaction kicks in, since s\ < 1/2. In this case, <S_ = S + = u~ x l 2 (\ — 2si)/2, and 
thus the 1.5PN coefficients in the relative equation of motion (1.4) take the form 

5 

A\.5PN — g<3, 

Bi. 5PN = ^Q, (7.9) 

where 

Q= T l ? (l-2 fll ) a =3^(l-2- 1 ) 9 . (7-10) 
At 2PN order, xa = $2 = 0, but S\ = Q ^ 0. In this case, the 2PN coefficients in (1.4) take the form 



. 9 Garni 
r 



a _ a gr , p GWi 

A 2PN — A 2PN ' *4 

r 

B2PN = B^pn ^ 2Q — — — . (7-11) 



Finally, the 2.5PN coefficients in Eq. (6.12) have the form 



01 = 3+^(9-277 + 3*/-), 

6 96 
25 

a 3 = -— Q(l-2?7 + V), 

h = l--j^Q{7-2r)-ty), 
96 

b 2 = 3-^Q(23-8t? + 3^), 

63 = ^Q( 13 + 2 V-W), (7.12) 
while the coefficients in the energy loss rate simplify to 

Ki - 12- yQ, 

k 2 = 11- jQ(17 + r?). (7.13) 

We find, somewhat surprisingly, that the motion of a mixed compact binary system through 2.5PN order differs from 
its general relativistic counterpart only by terms that depend on a single parameter Q, as defined by Eq. (7.10). 
Furthermore, all reference to the parameters Ai and A2, related to derivatives of the coupling function w(0), has 
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disappeared, in other words, the motion of mixed compact binary systems in general scalar-tensor theories through 
2.5PN order is formally identical to that in standard Brans-Dicke theory. The only way that a generalized scalar- 
tensor theory affects the motion differently than pure Brans-Dicke theory is through the value of the un-rescaled mass 
mi and the sensitivity s\ for a neutron star of a given central density and total number of baryons. 

The general conclusions reached in this paper about binary black holes and mixed binaries in scalar-tensor gravity 
were obtained from the near-zone gravitational fields. If these conclusions continue to hold for the gravitational- 
wave signal, then gravitational- wave observations of binary black holes will be unable to distinguish between general 
relativity and scalar-tensor theories, and observations of mixed black-hole neutron-star binaries will be essentially 
unable to distinguish between general scalar-tensor theories and Brans-Dicke theory. The radiative part of this 
problem, which will involve a derivation of the gravitational waveform to 2PN order, together with the energy flux, 
will be the subject of the second paper in this series. 
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Appendix A: Multipole moments for two-body systems 



Here we evaluate the multipole moments that appear in the radiation reaction expressions (6.6) and (6.7) to the 
order required to obtain 2.5PN-accurate contributions. The scalar dipole moment X\ in Eq. (6.6) must be evaluated 
to 1PN order. Substituting r s from Eq. (4.9e) and a s from Eq. (5.12c) to 1PN order into Eq. (3.7d), we obtain 



V s =GCmi4(l-2si) 



1 9 Gam 2 

2 1 r 
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(1^2). 



(Al) 



Most of the multipole moments that appear in the 2.5PN expressions (6.7) can be evaluated to the lowest PN order, 
so that we may write 



T 3 

j-ijk 

jgj 

•TO 
j-ijk 



G(l-C) ( 



m\x{ + m 2 x 2 



?) 



)• 



G(l - C) (m lX [ jk + m 2 xf 
G(l - C)e qab (rmvlxl 3 ' + m 2 v\x% 3 ' 
G( (mi(l - 2 Sl )x[ 3 + m 2 (l - 2s 2 )x\^ , 
G( (mr(l - 2 Sl )x[ jk + m 2 (l - 2 Sl )x$ k ) . 



(A2a) 
(A2b) 
(A2c) 
(A2d) 
(A2e) 



The exception to this rule is the scalar monopole moment M s = J M T s d 3 x; formally it contributes at 0.5PN order, as 
can be seen in Eq. (3.6d), but its leading contribution is constant in time, and hence it is the 1PN correction that 
matters. Inserting t s and a 8 from Eqs. (4.9e) and (5.12e) to 1PN order, we obtain 



M s = GCmi(l - 2si) 
Since the first term is constant, it can be dropped. 



Gam 2 



Pi 

1-4— 



+ (1^2). 



(A3) 
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